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ABSTRACT 

We derive the star formation histories of eight dwarf spheroidal (dSph) Milky Way satellite galaxies 
from their alpha element abundance patterns. Nearly 3000 stars from our previously published catalog 
(|Paper III ) comprise our data set. The average [a/Fe] ratios for all dSphs follow roughly the same 
path with increasing [Fe/H]. We do not observe the predicted knees in the [a/Fe] vs. [Fe/H] diagram, 
corresponding to the metallicity at which Type la supernovae begin to explode. Instead, we find that 
Type la supernova ejecta contribute to the abundances of all but the most metal-poor ([Fe/H] < —2.5) 
stars. We have also developed a chemical evolution model that tracks the star formation rate, Types II 
and la supernova explosions, and supernova feedback. Without metal enhancement in the supernova 
blowout, massive amounts of gas loss define the history of all dSphs except Fornax, the most luminous 
in our sample. All six of the best-fit model parameters correlate with dSph luminosity but not with 
velocity dispersion, half-light radius, or Galactocentric distance. 

Subject headings: galaxies: dwarf — galaxies: abundances — galaxies: evolution — Local Group 



1. INTRODUCTION 

Understanding the origins of galaxies requires under- 
standing the histories of their dark matter growth, gas 
flows, and star formation. Of these, the dark mat- 
ter growth is the most straightforward t o model (e.g., 
iDiemand et al.ll2007t [Springel et al.ll2008D . The gas flow 
history presents more difficult obstacles, such as colli- 
sional dissipation, gas cooling, stellar feedback, and con- 
version into stars. Despite the challenges, some models — 
built on top of dark matter simul ations — track all of thes e 
processes over cosmic time (e.g.. iGovernato et al.|[2007l ). 
The results of these models have observational conse- 
quences for the properties of the present stellar popula- 
tions of galaxies. 

1.1. Methods for Determining Star Formation Histories 

The star formation histories (SFHs) of galaxies may 
be deduced from the colors and magnitudes of the pop- 
ulation and from the spectra of the stars and gas, if 
present. Distant, unresolved galaxies display only a sin- 
gle, composite spectral energy distribution, which may 
be examined thr ough calibrations of spec trophotomet- 
ric indices (e.g., iGraves k, Schiavonl I2008D or, in some 
cases, spectral synt hesis dMcWillia m fe Bernsteinl 120081 : 
iColucci et al.1 12009). Nearer stellar systems may be re- 
solved both photometrically and spectroscopically. The 
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Hubble Space Telescope (HST) has enabled t he character- 
ization of the SFHs of many n e arby galaxies (iWeisz et al.l 
[20081 iDalcanton et~aTl [20091: iBernard et al.1 I2009D . in- 
cluding most of the dwarf galaxies in th e Local Group 
(jHoltzman et al.l [2001 lOrban et al.ll200l . 

Photometrically derived SFHs are most sensitive to 
young stars and metal-rich stars because the separation 
between isochrones increases with decreasing age and 
increasing metallicity. Elemental abundances obtained 
from spectroscopy do not give absolute ages, but they can 
provide finer relative time resolution f or old, metal-poor 
populations. IGilmore fc Wvsd (|1991f) showed that star 
formation bursts of varying duration and frequency in 
dwarf galaxies engrave signatures on the ratio of oxygen 
to iron as a function of metallicity. Because oxygen-rich 
Type II supernovae (SNe) explode within tens of Myr of a 
starburst, the oxygen content of stars forming soon after 
the burst will be high. Within hundreds of Myr, iron-rich 
Type la SNe begin to explode. The injection of iron into 
the interstellar medium (ISM) depresses the oxygen-to- 
iron ratio of subsequently forming stars. These processes 
are generalizable to other elements. The abundances of 
the next several elements with even atomic number be- 
yond oxygen — the alpha elements (Ne, Mg, Si, S, Ar, 
Ca, and Ti) — roughly scale with oxygen abundance. The 
abundances of iron-peak elements (V, Cr, Mn, Co, and 
Ni) roughly scale with iron abundance. The trend of 
the alpha-to-iron-peak ratio with iron-peak abundance, 
a proxy for elapsed time or integrated star formation, re- 
veals the relative star formation history with a resolution 
of about 10 Myr, the approximate timescale for a Type II 
SN. 

1.2. Chemical Evolution Models 

A glance at a diagram of [Mg/Fe] vs. [Fe/H] gives a 
qualitative sense of a galaxy's star formation history. 
Converting quantitative abundances into a quantitative 
SFH requires a chemical evolution model. iPagell (|1997l ) 
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described in detail how to create such a model, and 
iTolstov et al.l (|2009T ) reviewed recent progress on model- 
ing; th e SFHs of Local Group dwarf galaxies. iMatteuccll 
(2008) described the levels of approximation that the 
models assume. In general, more sophisticated and pre- 
sumably more accurate models reduce the number of ap- 
proximations. The most basic assumptions are instan- 
taneous recycling and instantaneous mixing. Considera- 
tion of stellar lifetimes and SN delay times removes the 
first approximation. Three-dimensional hydrodynamical 
simulations remove the second approximation. 

A chemical evolution model reflects the history not 
only of star formation but also of gas flow. A complete 
explanation of metallicity and alpha element distribu- 
tions requires both inflows and outflows. The metal- 
licity distribution functions (MDFs) of nearby Galactic 
G dwarfs cannot be explained with a closed box m odel 
(|van den Berghl [19621 ; ISchmidtl[l963l ). iPagel (fl997l) dis- 
cussed some of the proposed solutions to the G dwarf 
problem, including variable nucleosynthesis yields, bi- 
modal star formation, and pre-enrichment. On e of the 
most promising solutions is infalling matter (jLarsonl 
Il972[ ). Gases un doubtedly flow out of the galaxy, either 
from SN winds (jMathews fe Baked fl97lt iLarsonl [1971 
or s tripping from the influence of external or h ost galax- 
ies (jTinslev fe LarsonlfT979l: iLin fe Faberlfl98l . For ex- 
ample, interactions with the Milky Way could remove 
gas from the satellite galaxies discussed here. Both in- 
flows and outflows affect the star formation rate (SFR) 
throughout the history of the galaxy. Therefore, they 
shape the MDF and the trend of [a/Fe] with [Fe/H]. 

Chemical evolution models suffer from uncertain- 
ties in th e initial mass funct ion of stars and stellar 
lifetimes (IRomano et al.l 120051 ). nucleosynthesis yields 
(jRomano et al.1 120101 ), an d the delay time dist ribution 
(DTD) for Type la SNe ([Matteucci et al.ll2009| ). How- 
ever, these limitations have not prevented the mod- 
els from providing good fits to abundance data. Even 
models with some o f the first theoretical SN yields 
(jWooslev et al.l 119931 ) successfully reproduced the ob- 
served metalli city distribution and abund ance patterns 
in the Galaxy (Pagcl & Tautvaisienc 1995). Models with 
newer SN yields also match the sola r neighborhood abun - 
dance distributions very well (e.g.. IRomano et al.ll2010h . 
Nonetheless, uncertainties in the model assumptions do 
complicate the interpretation of the model results. For 
example, changing the Type la DTD, particularly the 
turn-on time, affects the derived timescale for star for- 
mation. The best way to circumvent these uncertainties 
is to apply the same model consistently to several systems 
and compare them differentially. Although the absolute 
ages or SFRs may be affected by systematic errors in the 
model, the relative quantities between different galaxies 
will be meaningful. 

Local Group dwarf galaxies make good subjects for 
chemical evolution models. First, the Local Group 
contains many reso lved dwarf galaxies (|Mateol 119981 : 
ITolstov et al.l 120091 ) with stars bright enough for 
medium- or high-resolution spectroscopy. Second, dwarf 
galaxies span a wide range of properties, including 
velocity dispersion and luminosity. The populations 
of the lowest luminosity galaxies enable the s t udy of 
star format i on on small scales ([Martin et al.l [2008a; 
iNorris et al.l [2008). The changes in populations for 



more luminous or more massive galaxi es show how 
star format i on res ponds to galaxy size (IMateol 1 1998; 
iKirbv et al.l 1201 0a'). Third, dwarf g alaxies host some 
of th e most metal - poor s tars known (IKirbv et al.1 | 200S 
2001 iGeha et al.l pol iCohen fc Huand 1 20091 I201C 
Frebel et all l2010allbl: iSimon et all 120101: INorris et alJ 
2010allbl : iStarkenburg et al.l 12010: iTafelmever et al.l 
l2010f ). and simlOb These stars retain the chemical 
imprint of the ISM when the Universe was less than 
1 Gyr old. Therefore, dwarf galaxies permit the study 
of star formation not only on small scales but also 
at early times. Finally, dwarf galaxies may be the 
primary building blocks for the Milky Wa y (MW) halo 
(|Searle fc Zirml Il978t IWhite fc Reesl [19781) . The stellar 
populations of the surviving dwarf galaxies may reflect 
the stellar populations of the dissolved building blocks, 
and they may show how the surviving satellites evolved 
since the time of rapid a ccretion onto the MW. 

In a series of articl es. lLanfranchi fc MatteuccH (|2003|. 
I2004J2007U2010D and lLanfranchi. Matteucci. fc Cescuttil 
(2006. 120081 ) presented numerical models that tracked 
the evolution of several elements in dSphs. The mod- 
els plausibly explained the MDFs and the available 
multi-element abundance measurements in dSphs. How- 
ever, large samples of published abundance measure- 
ments in any individual dSph have been spar s e until 
recently dShetrone et al.ll2009i: IKirbv et al.ll2009l l2010bl: 
ILetarte et al.l I2010D . Other chemical evolution mod- 
els of dSphs have e xamined the effects of reionization 
Fenner et al.1 I2006T) and sta r formation s tocha sticity 



Carigi fc Hernandez! l2008f ). iRecchi et al.1 ([2001) con- 



structed one of the first hydrodynamical models of dwarf 
galaxy evolution. In particular, they simulate d a galaxy 
similar to IZwl8. iMarcolini et al.1 (|2006l l2008| ) published 
hydrodynamical simulations of an isolated, Draco-like 
dSph. Their models relaxed the assumption of instan- 
taneous mixing and allowed inhomogeneous chemical en- 
richment. Some of the newest hydro dynamical models 
(|Revaz et al.1 120091: iSawala et aT1l2010D tracked both the 
kinematics and abundances of the stars as they form. 
They attempted to explain not only chemical abundance 
patterns but also dynamical properties of dSphs, such as 
the seemingly univ ersal dynamical mass measured w ithin 
their optical radii (|Mateolll998l; iStrigari et aT1l2008D and 
out t o the edge of their light distributions (jGilmore et all 
[200l . 

1.3. History of Chemical Analysis of Milky Way 
Satellites 

The earliest indications of heavy element abun- 
dance spreads among red giants of the dSph sys- 
tems in Draco, Ursa Minor, Sculptor, and Fornax 
were first obt ained b y the multi channel scanner ob- 
servations of iZinnl (19781. I1981D. i nitial efforts at 
spect r oscopy (|Norris fc Besselll 119781: IKinman fc Kraftl 
19801: IKinman et all 119811: IStetsonl 119841 : ISmithl 119841: 
Lch nert et al. I I1992D . and both broad and nar- 



row band photometry (jDemers. Kunkel. fc Hlrd^[l979l : 
iSmith fc Dopital 119831 ). The globular clusters of the 
Forn ax system proved to differ in their metallici- 
ties dZinn fc Perssonl 119811). The presence of carbon 
stars (lAaronson fc Mould! Il980t lAaronson et al.l 119821 
[19831: lAzzopardi et al.l Il985l)~ and so-called anomalous 
Cepheids (|Demarque fc Hirshfcld] 119751: INorris fc Zinnl 
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[T^IHirshfeldl[r980t[Snrrth fe Strvkeri 19861) further indi- 
cated the potential complexity of the stellar populations 
in dSphs. Carbon stars are exceedingly rare in globu- 
lar clusters, while the period-luminosity relations of the 
anomalous Cepheids implied t hat they are more m assive 
than typical cluster Cepheids (jZinn fe Searld lT976). As a 
consequence, by the mid-1980s, circumstantial evidence 
was building to suggest that dSphs had more complex 
and possibly more extensive star formation and chemical 
evolution histories than globular clusters. 

Since that time, the application of ground-based 
CCD and HST imaging has lead to greatly im- 
proved color-magnitude diagrams (CMDs) that have 
clearly shown the presence of significant internal age 
spreads within some of the Milky Way's retinue of 
dSphs . such as Car i na, F o rnax, Leo I, and Sextans 
/e.g.. iMighelll ITiM fl997t ISmecker-Hane et all ITiM 
Hurley- Keller. Mateo, fe Nemed 119981: [B uonann o et all 
1999t iGallart et al.l Il999allbl: ISaviane et al.l 120001: 
Lee et al.ll2009fl . Spectroscopy with large ground-based 
telescopes has demonstrated the presence of abun- 
dance inhomogeneities in t h e majority of these systems 
/e.g.. | Suntzeff et all 119931: ISmecker-Hane et al.l 1 1999 1: 
Shetrone et al.l I200U 120031 iTolstoy et all 120011 bood 

[Geisler et all 
2005aBbT 
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1.4. Chemical Evolution Models for the New Catalog 

In this article, we interpret the multi-element abun- 
dance distributions in eight dSphs with our own chemical 
evolution model. The data set is our catalog of abun- 
dances based on spectral synthesis of medium-resolution 
spectra fr om the DEIMOS sp ectrograph on the Keck II 
telescope (jKirbv et alll2010R Paper II) . The catalog con- 
tains 2961 stars with abundance measurements. The 
number of stars in each dSph ranges from 141 (Sextans) 
to 827 (Leo I). It is the largest homogeneous chemi- 
cal abundance data set in dwarf galaxies. The typical 
areal coverage is about 300 arcmin 2 at or near the cen- 
ter of each dSph. The median uncertainty on [Fe/H] is 
0.12 dex. The fraction of the sample with [Mg/Fe] un- 
certainties less than 0.2 (0.3) dex is 42% (53%). That 
fraction increases to 54% (69%) for [Ti/Fe], which is eas- 
ier to measure than [Mg/Fe]. For ([a/Fe]) (the average 
of [Mg/Fe], [Si/Fe], [Ca/Fe], and [Ti/Fe]), the fraction 
increases to 71% (88%). 

Our one-zone model is simple, but it incorporates some 
of the newest SN yields and the most recently measured 
DTD for Type la SNe. The biggest advantage of our data 
set is that it is homogeneous. All of the spectra were ob- 
tained with the same spectrograph configuration, and all 
of the abundances were measured with the same spectral 
synthesis code. Thus, the derived star formation and gas 
flow histories from our model — despite its simplicity — 
will be easy to interpret differentially. In other words, the 
absolute ages and star formation rates may be affected 
by model uncertainties, but the trends with galaxy prop- 
erties, such as luminosity, should reflect the true SFHs. 

We begin by describing our model (Sec. [5]). Then, we 
apply the model to the eight dSphs by finding the solu- 
tion that best matches the abundances. We discuss how 



our results compare to previous photometric and spec- 
troscopic studies (Sec. [3]). Next, we change some of the 
model variables to estimate the systematic errors in the 
derived SFHs (Sec. HJ). Then, we explore how the abun- 
dance distributions, SFHs, and gas flow histories change 
with galaxy properties such as luminosity and velocity 
dispersion (Sec. [5]). Finally, we enumerate our conclu- 
sions (Sec. [B]) . 

2. CHEMICAL EVOLUTION MODEL 

In order to pr ovide a ro ugh interpretation of the abun- 
dance trends in iPaper Ill 's catalog, we have developed a 
rudimentary model of chemical evolution. Table [TJ de- 
fines the symbol for each variable or constant in the 
model. The model supposes that a dwarf galaxy at any 
instant is a chemically homogeneous system that can 
accrete or lose gas. The ejecta of Type II SNe enrich 
the gas according to the total lifetime of massive (10 < 
M/M Q < 100) stars, while the T ype la SNe follow the 
observed DTD (M aoz et al.ll2010j. see belowl Stars form 
according to the iKroupa et alT " (fl993l ) initial mass func- 
tion (IMF, dN/dM = 0.31M" 2 - 2 for 0.5 < M/M Q < 1 
and dN/dM = 0.31A/~ 2 - 7 for M > 1 M Q ). 

The calculation tracks the mass of H, He, Mg, Si, Ca, 
Ti, and Fe at each time step (At — 1 Myr). The cal- 
culation is terminated when the system reaches zero gas 
mass. 

We define £,j(t) as the galaxy's gas mass of element j 
at time t. The galaxy's total gas mass at time t is 
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&e(t) H 
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20.4[£ Mg (i) 
f 1.07£ Fe (i) 



(2) 

The summation in Equation [1] is over all elements in the 
periodic table. However, our model tracks only seven el- 
ements. Therefore, we assume the ratio of the sum of 
all elements from Li to Ti, inclusive, to the sum of Mg, 
Si, Ca, and Ti is the same as i n the Sun. This ratio 
is 20.4 ([Anders fe Greve ssc 1989). Similarly, we assume 
the solar ratio for the sum of all elements V through 
Ge compared to Fe: 1.07. Elements beyond Ge are ne- 
glected. Equation [2] reflects these approximations. For 
convenience, we define the metallicity of the gas as fol- 
lows: 



Z = 



M gas (t) - - &e(«) 



(3) 



M gas (i) 

We also define the gas-phase mass fraction in an element 
3- 



(4) 



M sas (t) 

The following subsections explain the components of 
the models. Each component is expressed as the time 
change in £j(t), where £j = d£j(t)/dt. 

2.1. Star Formation Rate 

For simplicity, we assume that the star formation rate 
is a power law in the gas mass of the galaxy. With this 
assumption, 
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TABLE 1 

Chemical Evolution Model Variables 



Variable 


Description 


Units 


t 


Time since start of simulation 


Gyr 


M 


Mass of a single star 


Mq 




Gas mass in element j 


Mq 


Xj{t) 


Mass fraction in element j 


dimensionless 


Y 


Primordial helium mass fraction (Vjj c (0)) 


dimensionless 


M gas (t) 


Total gas mass 


Mq 


Z(t) 


Metal fraction (all elements heavier than He) 


dimensionless 


Ci(t) 


Time derivative of £j 


Mq Gyr- 1 


4>(*) 


Star formation rate, or rate of gas loss in element j due to star formation 


Mq Gyr" 1 


4-,n(t) 


Type II SN or HN yield rate for element j 


Mq Gyr" 1 




Fraction of HNe among stars with M > 20 Mq 


dimensionless 


O.nCM, Z) 


Mass of element j ejected by one Type II SN 


Mq 




Type la SN yield rate for element j 


Mq Gyr" 1 


^delay 


Type la SN delay time 


Gyr 


^la(tdelay) 


Type la SN delay time distribution 


SN Gyr" 1 Mq 


Cj.Ia 


Mass of element j ejected by one Type la SN 


Mq 


?j,AGbM 


AGB yield rate for element j 


Mq Gyr" 1 


Cj,AGB(M, Z) 


Mass of element j ejected by one AGB star 


Mq 


At, 


Normalization of star formation rate law (free parameter) 


Mq Gyr" 1 


a 


SFR exponent of Af ga s (free parameter) 


dimensionless 


A in 


Normalization of gas infall rate (free parameter) 


Mq Gyr" 1 


Tin 


Gas infall time constant (free parameter) 


Gyr 




Gas lost per SN (free parameter) 


Mq SN" 1 


Mgas(O) 


Initial gas mass (free parameter) 


Mq 



^' 3K> V.10 6 Mq 



(5) 



The variables A* and a are free parameters in the model. 
In the complete chemical evolution equation (Eq.[T6|). the 

sign of £j- » is negative because £j represents the gas mass, 
which is depleted due to star formation. 

E quation [5] is a generalization of a Kennicutt-Schmidt 
law ()Schmidti[l959[lKennicuttill998[ ). which connects the 
SFR to the gas surface density, £ gas - Surface density is 
perhaps more appropriate for disks than spheroids. De- 
siring a more three-dimensional property, we have used 
the gas mass, Af gas , instead of £ gas . The volume density, 
p gas , would be a better description, but the difference 
between Af gas and p gas is simply a constant because our 
model is one-zoned. 

2.2. Type II Supernovae 

In our model, stars more massive than 10 Mq and less 
massive t han 100 Mq explode according to their to tal 
lifetimes (|Padovani fc Matteuccilll993t lKodamalll997t ): 

r*(AT) = (l.2 (A//Af Q ) -1 - 85 + 0.003) Gyr (6) 

This formula is valid for stars more massive than 6.6 Mq, 
(inclusive of our entire m ass range for Type II SNe). 
iMaeder fc Mevnetl (|1989f ) give slightly different formu- 
las for stars less massive than 60 Mq, but the differences 
do not affect the chemical evolution model appreciably. 

Stars more massive than 100 Mq do not form in this 
model. The Type II SN ejecta are mixed homogeneously 
and instantaneously into the interstellar medium (ISM) 
of the entire dSph. 

We adopt the Ty pe II SN nucleosynthetic yields of 
iNomoto et all (|2006[ ). The symbol Q,ii(M, Z) represents 



the mass in element j ejected from the Type II SN explo- 
sion of a star with an initial mass M . It i s a function of 
both initial stellar mass and metallicity. INomoto et al.1 
tabulated the yields for seven initial masses ranging from 
13 Mq to 40 Mq and four metallicities from Z = to 
Z = 0.02. The total mass of the ejecta is always less 
than the birth mass of the star because the star loses 
some mass during its lifetime and because some mass is 
lo cked up foreve r in a SN remnant. 

INomoto et al.l modeled both normal core-collapse SNe 
and very energetic hypernovae (HNe). The lowest mass 
HN they model is 20 Mq. The fraction o f stars at least 
this massive that explode as HNe is £hn- INomoto et al.l 
adopted £h n = 0-5 for thei r own m odel of the solar neigh- 
borhood. iRomano et al.l (120 10; explored the cases of 
£hn = and 1. In our own experimentation, we have 
found that ehn = produces good matches to the dSph 
abundance patterns at the lowest values of [Fe/H], and 
we adopt this value for the model. In Sec. 14.21 we explore 
the effect of increasing 6hn on the model. 

The following integral gives the instantaneous change 
in gas mass from the ejecta of Type II SNe (Mq Gyr" 1 ): 

/■100 Mq 

4„=0.31M° 7 / Q >n {M,Z(t-n{M))) 



10 M 

x - TV (A/)) M- 2 - 7 dM 



(7) 



The coefficient 0.31 Mq 7 is the normalization from the 
IMF. This integral depends on the SN yields (£j,ii)) the 
recent star formation history (£*), and the high-mass 
IMF slope (Af- 2 7 ). In practice, this integral is per- 
formed numerically with Newton-Cotes integration over 
an array of 100 logarithmically spaced masses between 
10 Mq and 100 Mq. The values of £j,n and are inter- 
polated onto this array. The metallicity used to look up 
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Fig. 1. — Type la s upernova delay time distribution, as measured 
by Maoz ct al. ( 2010). The data come from a variety of star forma- 
tion environments, given in the figure legend. Equation[9]gives the 
expression for this function. Compare this figure to Mao z et alj s 
Fig. 2. 

the appropriate SN yields is consistent with the metal- 
licity of the gas at the time the exploding star formed. 
(In other words, at any given time step, the metallicities 
of the lower mass SNe are less than the metallicities of 
higher mass SNe from more recently formed stars.) 

The instantaneous Type II SN rate (SN Gyr -1 ) is given 
by a related integral: 



Nn = 0.31 M"- 7 



100 Mq 



t*(t-T*{M))M- A -' dM (8) 



10 M e 



This integral is performed over the same array of massive 
star lifetimes as a function of mass as for Eq.[7] The value 
will be used to determine the mass lost from SN winds 
(Sec.EU). 



2.3. Type la Supernovae 

We adopt the Type la SN yields of llwamoto et al.l 
(1999). The mass of element j ejected per Type la SN 
is C.j,ia- The SNe explode according to a function that 
approximat e s the delay time distribution observed by 
IMaoz et all (|2010t see Fig. [J). The following equation 
describes the adopted delay time distribution. 







(1 x 1CT 3 SN Gyr" 1 M" 1 ) 



Way < 0.1 Gyr 



/ tdclay ) 

V G y r ) 



idoiay > 0.1 Gyr" 



The variable ideiay is used instead of t to indicate that 
the DTD will be integrated from time t into the past. 

Unfortunately, the abundance distributions derived 
from the chemical evolution model depend sensitively 
on the normalization and turn-on time of $T a - Both 
of these quantities — particularly the turn-on time — have 
large uncertainties. The normalization affects [Fe/H] 
and the slope of [a/Fe] with [Fe/H]. We have chosen 
1 x 10~ 3 SN Gyr -1 M^ 1 for the normalization because 
that is the value that IMaoz et al.1 (|2010l ) reported. Even 



though the data in Fig. [T] are easily consistent with half 
that value, the larger value better reproduces the slope 
of [a/Fe] with [Fe/H] for many of the dSphs. The turn- 
on time determines the time or [Fe/H] at which [a/Fe] 
begins to drop. We have chosen 0.1 Gyr because that 
is approximately the maximum value acceptable for the 
DTD data (Fig. [[]). See Sec. 14.11 for a discussion of the 
effect of increasing this minimum delay time to 0.3 Gyr. 

The instantaneous Type la SN rate is given by com- 
bining ^i a with the past star formation history: 



£*(^dclay) ^la(t ~ ^dclay) dt delay 



(10) 



The mass returned to the ISM is the product of the SN la 
yields (Cj.ia) and the la rate: 



2.4. Asymptotic Giant Branch Stars 



(11) 



Winds from low- and intermediate-mass stars re- 
turn a small but significant amount of mass to the 
ISM. The stars lose less than 1% of this mass be- 
fore reaching the asymptotic gi ant branch (AGB, 
Ivan den Hoek fc Groenewegenl fT997) . Therefore, we con- 
sider mass loss on the AGB only. 

We adopt the AGB yields of iKarakad (|2010l ). who 
tracked all of the elements we consider here except Ca 
and Ti. (We assume that the fraction of Ca and Ti in 
AGB ejecta is the same as in the material that formed 
the star.) We assume all of the mass is ejected in the 
final time step of the star's lifetime. This assumption 
is appropriate because an AGB star's thermal pulsation 
period, during whi ch it loses most of its m ass, lasts on 
the order of 1 Myr (Marigo & Girardi 2007) , which is the 
length of one time step in our model. Equation [6] gives 
the lifetimes of sta rs more massive than 6 .6 Mq. Less 
massive st ars ob ey iPadovani k. Matteuccil 's (|1993l ) and 
iKodamaF s (|1997ft equation: 



0.334- 790 -0. 2232 [7. 764 -log ( i\//Af Q )] 

r»(M) = 10 5-TTT6 Gyr (12) 

Each AGB star ejects Ci.AGB solar masses of element j. 
Stars lighter than 10 Mq participate in AGB mass loss 
whereas stars heavier than 10 M© explode as Type II 
SNe (Sec. I2.2[) . The lower mass limit we consider for 
AGB stars is 0.865 Mq, which is the stellar lifetime cor- 
responding to the age of the Universe, 13.6 Gyr, accord- 
ing to Eq. [12j The AGB mass return rate in Mq Gyr -1 
is given by 



1 Mq 

0.865 Mq 



4aGB = 0.31M£ 2 / 6,AGB(M,Z(t-T„(M))) 



x f„(t - r*(M)) AT 



dM 



0.31 M, 



© 



10 Mq 



1 Mq 



Q, AGB (M,Z(t-n(M))) 



x £*(t-T*(M))M~ 



dM 



(13) 



Compared to SN ejecta, AGB ejecta affect the chemical 
evolution of the elements considered here to a small de- 
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gree. AGB ejecta are more important for other elements, 
such as C, N, and O. 

2.5. Gas Infall 

Infall of gas during the star form ation lifetime of a 
dSph is required to explain its MDF (jKirbv et al.ll2010al 
Paper III). Therefore, our model allows pristine gas to 
fall into the dSph. The gas has a helium fraction of 
Y = Xho(O) = 0.2486, which is t he value obtained when 
the WMAP7 ((Larson et al.ll2010D barvon-to-p hoton ratio 
is applied to the formula of iSteigmanl ([2007). The rest 
of the infalling gas is hydrogen. 

The MDFs of the dSphs are generally more peaked 
than a closed box model predicts. One scenario that 
explain such a distri bution is gas infall that first incr eases 
and then decreases (|Lvnden-Belllll975t lPagellll997t ). We 
find that a quick increase of the rate of gas falling into 
the galaxy followed by a slower decrease in the infall rate 
does well at reproducing the data. We parametrize the 
gas infall rate as follows. 

£ jM = A Sa X j (t = 0) e-^" (14) 

The term Xj(t — 0) means that the infalling gas is pri- 
mordial (metal- free) . The variables A ln and t ivl are free 
parameters in the model. 

2.6. Supernova Winds 

The MDFs of dSphs require gas outflow. If that were 
not the case, the metallicities would approach the super- 
nova yields, which are much larger than observed in even 
the most metal-rich star in any dSph. Gas may be lost 
through supernova winds, stellar winds, or gas stripping 
from an external source. All of these sources undoubt- 
edly occur over a dSph's lifetime, but supernova winds 
are the most straightforward to include in a chemical 
evolution model. We ignore other sources of gas loss. 

Our computation of gas loss is fairly simple. The 
galaxy loses a fixed amount of gas for every super- 
nova that explodes. The blown-out gas mass does not 
vary with SN type because the e xplosion energies fo r 
Types II and la SNe are si milar. SeelRecchi et al.l (12001'). 
iRomano et ail (|2006h . and iMarcolini et al.l (|2008T ) for ex- 
amples of chemical evolution models that treated the en- 
ergy input from the two SNe types differently. The rate 
of gas loss is 

4out = A out Xj (N u + Nied (15) 

The parameter A ont is a free parameter in the model. An 
energy argument shows that the ejected gas mass is of 
the order of 10 4 Mq SN -1 . O ne supernova explodes wit h 
a typical energy of 10 51 erg (|Wooslev fc Weaver! 119951) . 
In the late stages of expansion , the kinetic energy of the 
ejecta is E ci ~ 8.5 x 10 49 erg ((Thornton et al.1 1 1998( 1. A 
typical line-of-sight velocity dispersion for a dwarf galaxy 
is <7i os ~ 10 km s . Given the virial theorem (GM/R = 
3of os ) and the escape velocity (v^ sc — 2GM/R), then 
the gas mass ejected as a result of SN blowout is M e j = 
E ci /v 2 csc = E ej /(6afJ ~ 7 x 10 3 M SN" 1 . 

A metal-enhanced wind can prevent the galaxy from 
becoming too metal-rich without such a large gas loss 



(|Vaderlll986l) . For simplicity, we assume that the SN 
winds have the same chemical content as the gas remain- 
ing in the galaxy. See Sec. 14.31 for a further discussion of 
including metal-enhanced winds in the model. 

2.7. Complete Chemical Evolution Equation 

The complete equation that describes the chemical evo- 
lution of the galaxy's gas is 

0(t) = M gas (0)+ f (-4* + 4mi + 4ia + (16) 
Jo 

£j,AGB + £j,in — £i,out) dt 

The initial gas mass, M gas (0), is a free parameter. A 
non-zero initial gas mass may seem inconsistent with 
Eq.[5]because the gas should form stars as it falls into the 
galaxy. However, the galaxy could acquire gas available 
for star formation — via gravitational collapse or cooling, 
for example — on a timescale faster than the star forma- 
tion timescale. We will show that the non-zero initial gas 
mass is more important for the more luminous dSphs. 

2.8. Shortcomings of the Model 

Our model incorporates realistic conditions in dwarf 
galaxies. We mode l chemical evolutio n using an observed 
Type la SN DTD (jMaoz et al J 12010( 1 . We also take into 
account the lifetimes of Type II SN progenitors, rather 
than assuming instantaneous recycling. The delay helps 
to shape the metal-poor abundance distributions because 
it affects the rapid rise in metallicity after the onset of 
star formation. 

However, our model is not as sophisticated as some 
other chemical e volut i on models of dwarf galaxies (e.g., 
Marcolini et~aTl 120081: iRevaz et all 120091: iSawala et all 
20101) . In the next section, we show the best model fits 
to eight different MW satellite galaxies. The simplicity 
of our model reduces the computational demand of find- 
ing the best solution. Nonetheless, we enumerate some 
shortcomings which affect the interpretation of the abun- 
dance distributions. 

1. The turn- on time for Type la SNe is poorly con- 
strained. iMaoz et "all (|2010D showed that it is al- 
most certainly 0.1 Gyr or less (at least in the Mag- 
ellanic Clouds and higher redshift elliptical galax- 
ies), but the DTD slope (t<jeiay) ^ s divergent as 
£deiay approaches zero. Therefore, the number of 
Type la SN that explode shortly after their progen- 
itors form depends sensitively on the turn-on time. 
The uncertainty in the turn-on time translates to 
a large uncertainty in the Fe abundance distribu- 
tion. With all other model parameters held fixed, 
an earlier turn-on time would cause the metallicity 
of the MDF peak to increa se and [a/Fel at a given 
metallicity to decrease. SeelMattcu cci et al.l {2009) 
for a detailed discussion of the effect of adjusting 
the ratio of prompt to delayed Type la SNe. 

2. The SN yields are imperfect. As we mention in 
Sec. [3J we needed to increa se the [Mg/H] outpu t 
of the model by 0.2 dex (see iFrancois et al"1 l2004'). 
Furthermore, Ti is severely underproduced in our 
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model. Therefore, we do not consider Ti abun- 
dances at all. 

3. Our model assumes instantaneous mixing. Re- 
laxing this approximation would require multiple 
zones, which we do not consi der for the sake of 
computational s implic i ty. S e e iMori et al. 



(2002) 



iMarcolini et~aH (12 006. 200|), iRevaz et al.l ([2009T ) 
and iSawala et al.l (|2010f ) for three-dimensional 
chemical models of dwarf galaxies. 

4. We also assume instantaneous gas cooling. The 
cooling time for gas to become available for star 
formation (after accretion or ejection from SNe and 
AGB stars) may be longer than the model time 
step, At = 1 Myr. A more proper treatment of the 
cooling time, such as in a hydrodynamical model, 
might result in slightly longer SF durations that we 
derive with instantaneous cooling. 

5. On a related note, we also ignore dynamical pro- 
cesses. Our adoption of a single value of A out , the 
gas ejected from the galaxy in the wind of one su- 
pernova, implicitly assumes that the potential of 
the galaxy is homogeneous and static. This as- 
sumption is inconsistent with our allowance of gas 
to flow into the galaxy. Although dark matter dom- 
inates the dynamical mass of dSphs, they undoubt- 
edly change their dark mat ter masses durin g their 
star formation lifetimes dRobertson et al.l 120051 : 
iBullock fc Johnston! [20051 : Uohnston et al.1 120081 ) . 
Furthermore, baryonic (adiabatic) contraction can 
affect star formati on and feedback in the d ense cen- 
ters of the dSphs ( N apolitano et al]|2010f ). 

6. We consider only one parametrization of the gas in- 
fall rate. Because the star formation rate is propor- 
tional to the gas mass, the gas infall rate essentially 
shapes the differential MDF. Differently shaped gas 
infall histories might better reproduce the dSph 
MDFs. External influences on the gas flow (or al- 
ternatively, availability of gas cool enough to form 
stars) that we do no t consider include reionization 
(IBullock et all I2000D a nd tidal and ram pressure 
stripping (|Lin fc Faber|[l983fl . 



7. We model only one episode of star formation. 
CMDs have revealed extended and possibly bursty 
SFHs in several dSphs in our sample (Fornax and 
Leo I and II). These bursts will not be included in 
our model. In these cases, we defer to the photo- 
metrically derived SFHs. In fact, we suggest for fu- 
ture study a more sophisticated analysis that mod- 
els both the CMD and abundance distributions. 

8. The infalling gas is assumed to be metal-free at 
all times. In reality, the metallicity may have in- 
creased over time because the source of the new 
gas may have been blowout from prior SF episodes 
in the galaxy in question or other galaxies. This 
gas would have been enriched by SNe and other 
nucleosynthetic sources. 

9. The modeling result for a given galaxy repre- 
sents only part of that galaxy's stellar popula- 
tion. Our spectroscopic samples were centrally 



concentrated to maximize the number of mem- 
ber stars on a DEIMOS slitmask, but most dSphs 
have radial populati on gradients (e.g., Sculptor, 
Battaglia et al.l 120081 ). As a result, we preferen- 
tially probe the younger, more metal-rich popu- 
lations. MW satellite galaxies al so shed stars as 
they in teract with the Galaxy. Maicwsk i et al.l 
(2000b) identified stars fro m the Carina dSph be- 
yon d Carina's tidal rad ius. Maicws ki et al.l (2002) 
and lMufioz et al.l (2006) discussed the implications 
for Carina's present stellar population. In particu- 
lar, the remaining stars are on average younger and 
more metal-rich than the lost stars. Consequently, 
the spectroscopic sample favors the younger, more 
metal-rich stars. 

Some of these shortcomings are observational or the- 
oretical uncertainties (1-2), which can only be resolved 
with a more thorough investigation of SN rates or yields. 
Others are simplifications (3-8), which can be resolved 
with more sophisticated models. The last shortcoming 
(9) could be resolved by an intensive, wide-field campaign 
with the intent to recover spectra for a magnitude-limited 
sample of red giants in a dSph. This project would re- 
quire a great deal of telescope time, but it could be ac- 
complished in principle for one or two dSphs. Foreground 
contamination could be minimized by selecting a dSph 
at high Galactic latit ude or photometrically pre-selecting 
likely members (e.g.. iMaiewski et al.ll2000al ). 

3. GAS FLOW AND STAR FORMATION HISTORIES 

We apply our chemical evolution model to eight dSphs: 
Fornax, Leo I, Sculptor, Leo II, Sextans, Draco, Canes 
Venatici I, and Ursa Mino r. We use the abundance mea- 
surements from iPaper iB . For each galaxy, we attempt 
to match simultaneously the distribution of [Fe/H] and 
the trends of [Mg/Fe], [Si/Fe], and [Ca/Fe] with [Fe/H] 
by adjusting the six free parameters listed at the bottom 
of Table [2 

Unfortunately, some elemental abundances could not 
be matched for any combination of parameter values. In 
particular, the model underpredicts [Mg/H] and [Ti/H]. 
iFrancois et al.l ([2004D constructed a chemical evolution 
model for the Milky Way and also encountered trouble 
in reproducing the yields. They concluded that the SN 
yields should be modified. They specifically singled out 
Mg for being underproduced by both Type la SNe and 
low-mass Type II SNe. We feel comfortable modifying 
the model results for [Mg/H] because chemical evolu- 
tion models by different authors over a wide range of 
galaxy masses and ages indicate that such modification 
is necessary. We add 0.2 dex to [Mg/H] to bring the 
mo del into better agreem ent with the data. However, 
the iNomoto et all (|2006h Type II SN yield for [Ti/Fe] 
is about —0.1 dex, which is far below the value observed 
for metal-poor stars in dSphs or in the MW halo. Rather 
than attempting to corre ct such a large deficit, we ignore 
the model result for Ti. INomoto et al.l also ignore their 
Ti yields in their own chemical evolution model of the 
sola r neighborh ood. 

In IPaper IIP , we found the best- fit analytical chemical 
evolution models for the same eight dSphs based on their 
MDFs alone. We repeat t he process here for our more 
sophisticated model. As in IPaper III we use maximum 
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likelihood estimation to find the best-fit model parame- 
ters. 

The likelihood that a particular model matches the 
data is the product of probability distributions. Each 
star is represented by a probability distribution in a 
four-dimensional space. The four dimensions are [Fe/H], 
[Mg/Fe], [Si/Fe], and [Ca/Fe]. We denote these quanti- 
ties as £ij, where i represents the i star and j identifies 
one of the four element ratios. The Gaussian is centered 
on the star's observed values. The width in each axis is 
the estimate of measurement uncertainty {Seij) in that 
quantity. Stars with larger uncertainties have less weight 
in the likelihood calculation than stars with smaller un- 
certainties. (Although Figs. 2-9 show only stars with 
uncertainties less than 0.3 dex, there is no error cut in 
the likelihood calculation. Instead, we downweight stars 
with large uncertainties.) The chemical evolution model 
traces a path ej(t) in the four-dimensional space. The 
probability that a star formed at a point t on the path is 
dP/dt = M*(t)/M*, where M* is the galaxy's final stel- 
lar mass. The likelihood that one star conforms to the 
model is the line integral of dP/dt along the path ej (t) . 
The total likelihood L is the product of the individual 
likelihoods of the N stars: 



1 



— (Af*,obs — M* iInoc iel) 2 



dt 



'27n5M* i0 bs 
f 



exp • 



- ( iw gas obs 

2-7T 6 Af gaSi0 bs 2(dM gaSi0bs 



2(<5M*, obs ) 2 

2 \ 0.1N 



(17) 



The second line of the equation requires that the final 
stellar mass of the model (M* jmo dci) matches the ob- 
served stellar mass (M* j0bs ) within the observ ational un- 
certai nties. We adopt the stellar masses of Woo et al.l 
(2008). They did not study Canes Venatici I. We as- 
sume that galaxy has about the same stellar mass as Ursa 
Minor because it has the same luminosity within the ob- 
servational uncertainties. The third line of the equation 
assures that the dSph ends up gas free. We fairly ar- 
bitrarily assume an uncertainty of #M gaS)0bs = 10 3 M 
because even lower values of <5M gaSj0bs cause the chem- 
ical evolution model to converge on spurious solutions. 
The exponent 0.1 N sets the relative influence of the fi- 
nal stellar and gas mass compared to the abundance dis- 
tributions. This value was chosen so that these quanti- 
ties did not dominate the likelihood but also so that the 
modeled galaxies ended up gas-free and with about the 
correct stellar mass. 

For computational simplicity, we minimize the quan- 
tity L = — IriL: 



(-^gas,obs)' ! 

2{8M„ 



i.obs ) 



+ ln(27r) + ln(5M*, obs ) + ln(<5M gaS!0bs ( 



We find the values of the six parameters that minimize L 
using Powell's method. We calculate uncertainties on the 
model parameters via a Monte Carlo Markov chain. We 
perform at least 10 4 trials for each dSph after a burn-in 
period of 10 3 trials. The dSphs with shorter SF dura- 
tions require less computation time, and we were able to 
pe rform up to 5 x 10 4 trials for some of the dSphs. As 
in iPaper III the model uncertainties are the two-sided 
68.3% confidence intervals. These uncertainties incorpo- 
rate only observational uncertainty and not systematic 
model errors. Table [5] lists the solutions for each dSph in 
order of decreasing luminosity. 

Table [3] lists the total star formation durations for the 
most likely models. The duration is not a free parame- 
ter but a result of the model. The table also lists some 
timescales derived from HST CMDs (|Dolphin et all20"05t 
IQrban et al.l 12008). It is not possible to measure pho- 
tometrically the total star formation duration for pre- 
dominantly ancient stellar populations because 10 Gyr 
isochrones are extremely similar to 13 Gyr isochrones. 
Therefore, we have quoted /iog, the fraction of stars 
formed more recently than 10 Gyr. For small or zero 
values of fioc the CMD shows that the population is an- 
cient, but there is no time resolutio n. We also show th e 
stellar mass-weighted mean age r (|Orban et al.l [20081. 
For the three dSphs with intermediate- aged populations 
(Fornax and Leo I and II), r combined with fxoG gives 
some idea of the star formation duration. For example, 
Fornax formed 1 — /iog = 27% of its stars beyond 10 Gyr 
ago, but the mean age is just 7.4 Gyr. Half of Fornax's 
stars formed over at least 2.6 Gyr, and the other half 
formed even more recently. Our abundance-derived du- 
ration of 1.3 Gyr is inconsistent with this photometric 
star formation duration. For Fornax and Leo I and II, 
we defer to the photometrically derived SFHs (see item 7 
of Sec. 12.81) . They are more realistic because they per- 
mit an arbitrary number of SF episodes. For the galaxies 
whose CMDs identify them to be ancient, our abundance 
distributions are far more sensitive probes of the SF du- 
ration than the CMD. 

In the following sections, we discuss the derived star 
formation and gas flow histories for each dSph and com- 
pare them to previous photometrically and spectroscop- 
ically derived SFHs. 

3.1. Fornax 



N 

£ 

z=l 



111 



0.17V 



n 



exp 



DdclJ 



2(£A/*,obs) 2 



We begin our discussion with the most luminous of the 
mostly intact MW dSph satellites, Fornax. Its [a/Fe] dis- 
tribution (Fig. [5]) shows the least evidence of correlation 
with [Fe/H] of all eight dSphs studied here. In the range 
— 1.3 < [Fe/H] < —0.5, the four [a/Fe] element ratios 
span almost 1 dex at a fixed metallicity with no evidence 
of a slope with [Fe/H]. The rarer stars more metal poor 
than [Fe/H] w — 1.3 have higher average [a/Fe]. 

The large range of [a/Fe] and the lack of correlation 
with [Fe/H] each suggest bursty or inhomogeneous star 
— (ejj — £j(i)) 2 \ M*(t) formation. A bursty SFH would cause spikes and depres- 
sions in [a/Fe] as [Fe/ H] increases monotonically (e.g., 
IGilmore fc Wvsl 11991 ) . even if the star formation were 
well-mixed over the whole galaxy at all times. Measure- 
ment uncertainties might blur the division between the 
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TABLE 2 

Galaxy Properties and Chemical Evolution Model Parameters 



dSph L M, <ti os R e Dqc A* 

(10 5 Lq) (1O«M ) (km 8-1) ( P c) (kpc) (%^) f 15 ^) ( G ^) I — 



Fornax 


180 ± 


50 


190 ± 50 


10.7 


± 


0.2 


714 


± 


40 


141 


± 


12 


Leo I 


56 ± 


16 


45 ± 13 


9.0 


± 


0.4 


295 


± 


49 


257 


± 


76 


Sculptor 


22 ± 


10 


12 ±5 


9.0 


± 


0.2 


282 


± 


41 


85 


± 


23 


Leo II 


6.6 ± 


1.9 


14 ±4 


6.6 


± 


0.5 


177 


± 


13 


221 


± 


50 


Sextans 


4.1 ± 


1.2 


8.5 ±2.4 


7.1 


± 


0.3 


768 


± 


47 


98 


± 


13 


Draco 


2.7 ± 


0.4 


9.1 ± 1.4 


10.1 


± 


0.5 


220 


± 


11 


92 


± 


29 


Can. Ven. I 


2.3 ± 


0.4 


6±2 


7.6 


± 


0.5 


546 


± 


36 


210 


± 


29 


Ursa Minor 


2.2 ± 


0.7 


5.6 ± 1.7 


11.5 


± 


0.6 


445 


± 


14 


70 


± 


19 



5 02+ 2 ' 21 98+ - 15 2 46+ ' 70 31+ - 01 1 5 i+0 03 

°- uz -1.00 u - a5 -0.04 Z -* D -0.17 "■ ol -0.04 1 - ol -0.06 

n Q9 +0.67 n71 +0.01 i 17 +0.05 n oc+0.02 o on+0.16 

u - 3 -0.05 u -' 1 -0.17 1 - 1 '-0.10 u - od -0.01 °- o; '-0.08 

0.47±°;» 0.83^ 0.70±g;Sl 0.27t° f 2 5.361°;!? 0.50^ 

0.431H? Q.eA 0.48±°$ 0.42«f 9 6.59±°;g 0.05^ 

0.52l°;« 0.50l°-j° 1.15±g;S 0.22«°l 9.60t°;?° 1.55j 

0.881^0 0.34l°;« 1.27±g : ?| 0.221°;°! 9.5ll°;« 2.32^ 

0.46l°;« 0.361°;^ 0.86l°;g 0.2ll°;°* 8.83l°;f° 0.27^ 

l-21±8:ff 0.26l°;° 7 1.47lg;?| 0.17l°;° 2 11.04±g;S 0.54^ 



References. — L (luminosity): Martin ct al. (2008a) for Canes Vcnatici I, Irwin & Hatzidimitriou (1995) for the others. M* (stellar mass): Woo ct al. (2008), 
except t hat we have assum ed that Canes Venatici I has about the same M* as Ursa Minor. cti ob (line- of-sight veloci ty dis persion) a nd R e (2-D projected half-light 
radius): Wolf ct al. (2010) and references therein. Dqc (Galactoccntric distance): Coordinates from Mateo ( 1998). Sec Paper II, Table 1, for the sources of the 
heliocentric distances. 



TABLE 3 
Star Formation Durations 



dSph 


Duration 3 
(Gyr) 


f 10G (D05) b 


hoc (O08) c 


r (O08) 
(Gyr) 


Fornax 


1.3 


0.73 


0.73 


7.4 


Leo I 


1.4 


0.75 


0.76 


6.4 


Sculptor 


1.1 


0.05 


0.14 


12.6 


Leo II 


1.6 


0.56 


0.70 


8.8 


Sextans 


0.8 


0.00 


0.00 c 


12.0 


Draco 


0.7 


0.06 


0.49 


10.9 


Can. Ven. I 


0.9 








Ursa Minor 


0.4 


0.00 


0.00 e 


12.0 



Note. — Our star formation durations for Fornax and Leo I and II arc 
almost certainly too short because our chemical evolution model does not 
permit multiple SF bursts. 

a Star formation duration derived from our model, based on spectroscopic, 
multi-clement abundances. 

k Fraction of stars formed more recently than 10 Gyr ago, based on an 
analysis of HST photometry JDolphin e"t~aT||2005f) . 

c Fraction of stars formed more recently than 10 Gyr ago, based on a 
different analysis of HST photometry (Orban et alj|25oa >. 
^ Stellar mass-weighted mean age, based on Hubble Space Telescope pho- 
t ometry llOrban et alll2008l) . 

e IQrban et~aT1 (2008) did not measure these values but took them from 
IDolphin et al.l 120051) . 



[a/Fe] spikes in different bursts. Alternatively, if the 
SN nucleosynthetic products were not well-mixed, the 
[a/Fe] value of a star would reflect the particular SFH of 
its birth site rather than the galaxy as a whole. Conse- 
quently, the abundance distribution would be a compos- 
ite of several different SFHs. Coupled with measurement 
uncertainties, the composite distribution may look like 
an uncorrelated scatter of points, such as the distribu- 
tion in Fig. [2] Burstiness and inhomogeneity are not 
mutually exclusive. Both processes might have affected 
Fornax's SFH. 

Based on HST/ Wide Field Planetar y Cam era 2 
(WFPC2) photometry, IBuonanno et al.l (|1999f ) sur- 
mised that the field (not globular cluster) population 
of Fornax endured three major bursts of star for- 

iSaviane et al.l 



matio n separated by about 3 Gyr. 
(2000 ). [Battaglia et all (1200(1 . iGullieuszik et al.l (|2007|). 
and I Coleman fc de Jond (|2008T l provided additional 
photometric and spectroscopic evidence of multi- 



ple di screte population s, incl ud ing a burst 4 Gyr 
ago. IGrebel fc Stetson! ()1999D . IBattaglia et al.l . and 
I Coleman fc de Jongl additionally showed that the 
younger, more metal-rich populations are more centrally 
concentrated. Thus, it seems that star formation in 
Fornax was both bursty and inhomogeneous. 

Our chemical evolution model is incompatible with 
Fornax's complex SFH. First, we model the SFR as a 
smooth function, not a bursty one. Second, the model 
has only one zone and does not account for spatially 
segregated star formation. Consequently, the SFH de- 
rived from our model should be viewed with skepticism. 
Most notably, we derive a total star formation duration 
of 1.3 Gyr (the time at which star formation and SN 
winds exhausted the gas supply, thereby truncating star 
formation), whereas every photometric study shows that 
star formation in Fornax lasted for most of the age of 
the Universe. In addition, the model does not match the 
observed flatness of the [a/Fe] distribution for the bulk 
of the stars. However, the model does share one im- 
portant quality with photometrically derived SFHs: The 
initial metal enrichment is very rapid. The metallicity 
in our model reaches [Fe/H] = — 1 at 0.3 Gy r after the 
commencement of star formation. iPont e t al. (2004|) de- 
duced that Fornax reached [Fe/H] = — 1 within a few 
Gyr. One advantage of a spectroscopically derived SFH 
is that it is sensitive to relative ages, whereas a photo- 
metrically derived SFH is sensitive to absolute ages but 
h as poor age re s olution for old populations. 

iLetarte et alj (|2010f ) measured multi-element abun- 
dances from higher r esolution spectra of 81 Fornax mem- 
bers. We showed in IPaper ij that our abundance mea- 
surements match theirs very well. They pointed out 
that centrally selected stars in Fornax will preferentially 
sample the young, metal-rich component. In fact, the 
most metal- poor star k nown in Fornax ([Fe/H] = —3.66, 
lTafelmeveret~alT2 010) is very far (43') from the center of 
the dSph. The discovery emphasizes that selecting stars 
in the center of the dSph biases the age and metallicity 
distribution. 

3.2. Leo I 

Leo I is the second most massive dwarf galaxy in our 
sample. The [a/Fe] distribution of Leo I shows a mod- 



10 



Kirby et al. 



EEs. 

T3 



E5j 

so 



'ir, 



U 





t (Gyr) 



Fig. 2. — The observed abundance ratios and the best-fit gas flow and star formation history model for Fornax. Left: The top panel 
shows the observed MDF as the black histogram and the modeled MDF in red. The model is convolved with an uncertainty function to 
mimic the broadening of the histogram induced by observational error. A cross marks each 1 Myr time step, but these are too closely 
spaced to discern for most of the metallicity range. Very few stars are expected to have formed at the low metallicities where the crosses are 
distinguishable. The other panels show the observed [Mg/Fe], [Si/Fe], [Ca/Fe], and [Ti/Fe] ratios as black points whose sizes are inversely 
proportional to measurement uncertainties. Only points with uncertainties less than 0.3 dex are shown. The red lines show the abundance 
ratios of the stars and gas at each time step. We do not show the model results for [Ti/Fc] because the SN yields are inaccurate. Right: The 
gas flow and star formation history for the best-fit model. From top to bottom, the panels show the gas inflow rate; the stellar, gas-phase, 
and total baryonic mass; the star formation rate; the iron and magnesium abundances; and the [Mg/Fe] ratio, all as a function of time. 
The second panel also gives the final stellar mass in the model. 



erate correlation with 
metallicity stars ([Fe/H 



Fe/H]. In particular, the lower 
< —1.5) show on average higher 
[ q/Fel (except for Ti) than the more metal-rich stars. 

ILee et alj (fl993h obtained the first CCD-based CMD 
of Leo I, and they found hin ts of a young (3 Gyr) pop- 
ulation. iCaputo et all (|1999Q and iGallart et all (|1999ah 
conducted the first comprehensive studies of Leo Fs SFH 
using CMDs obtained with HST/WFPC2. Because these 
CMDs reached the main-sequence turnoff of the oldest 
(> 10 Gyr) populations, they were able to study the 
multiple stellar populations and complex SFH. Leo I was 
thought to be unique among the MW satellite dSphs 
for lacking a conspicuous horizontal branch (HB) until 
a 12' x 12' ground-based survey of Leo I by I Held et ahl 
(2000) revealed a HB structure in its CMD. The exis- 
tence of both an extended blue HB and RR Lyrae stars 
(|Held et alj |2001[ ) suggested that Leo I is in fact simi- 
lar to other local dSph galaxies in having a > 10 Gyr 
population, but the majority of stars were still believed 
to have formed later than 7 Gyr ago. However, a re- 
cent CMD obtained with H ST/ Advanced Camera for 
Surveys/ Wide Field Camera (|Smecker-Hane et aLll2009| ) 
reached far deeper than the earlier ones and showed that 
at least half of the stars were in fact formed more than 



9 Gyr ago, which is c onsistent w i th the abundant RR 
Lyrae stars found b y iHeld et ahl ((2001). In addition, 
iSmecker-Hane et al. | combined t heir C MD with the spec- 
troscopic MDF of iBosler et all (|2007j ) to find that Leo I 
experienced two episodes of star formation around 2 and 
5 Gyr ago. 

Because our chemical evolution models halt when the 
gas mass drops to zero, we are unable to recover the 
later phases of SFH (i.e., the two bursts at 2 and 5 Gyr 
ago). Nonetheless, our model provides insights into 
the early phase with better time resolution. Overall, 
our model matches the observed trend of [a/Fe] with 
[Fe/H] fairly well, but the model MDF slightly overpre- 
dicts the frequency of metal-rich stars. The observed 
MDF also shows a more pronounced peak at [Fe/H] 
= —1.4 than the model. The initial starburst that 
likely led to the formation of Leo I lasted for about 
1.4 Gyr. This is much shorter than the star forma- 
tion duration of ~ 5 Gyr derived by photometric stud- 
ies. As with other galaxies in our sample, adding bursti- 
ness to our model would help re solve these discrepancies. 
lLanfranchi &; Matteuccil (|2010f ) suggested that Leo I is 
characterized by a low SFR and intense galactic wind. 
The main difference between their model and ours is 
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Fig. 3. — The observed abundance ratios and the best-fit gas flow 
explanation. 

that we start with a much higher gas mass (by a fac- 
tor of ~ 400). Also, our model requires a highly efficient 
SFR to match the observ ed MDF. The discrepancies with 
ILanfranchi fc Matteuccil partly result from our choice to 
use unenhanced galactic winds. Metal-enhanced winds 
would reduce the amount of gas required to be blown 
out. As for Fornax, our model is qualitatively consis- 
tent with previously derived SFHs in the sense that the 
overall metallicity increases quickly at early times. 

Leo Fs orbital dynamics, a s studied by ISohn et aLl 
(|2007fl and iMateo et all (|2008| ). indicate close passes to 
the center of the MW. The dSph almost certainly lost 
stars in tidal interactions near its perigalacticon. The 
prevalence of an intermediate- aged (rather than old) pop- 
ulation in Leo I may be a consequence of this tidal strip- 
ping. Because the stripped stars do not fall in our spec- 
troscopic sample, our model does not represent some 
stars that formed early in Leo Fs history (see Sec. 12.81 
item 9). 

3.3. Sculptor 

Our chemical evolution model for Sculptor produces 
one of the best fits to the abundance distributions (Fig.[4j) 
out of all of the dSp hs, particularly for the asymmetrical 
MDF. In iPaper III we could not reproduce the width 
of Sculptor's MDF with an analytical model of chemical 
evolution. Our more sophisticated model, which more 
properly treats Fe as a secondary nucleosynthetic prod- 
uct with multiple origins (Types II and la SNe) , yields a 
broad, well-matched MDF for the appropriate choice of 




t (Gyr) 

and star formation history model for Leo I. See Fig. [2] for a detailed 

parameters. The combination of a low SFR normaliza- 
tion (A*) and low initial gas mass maintains a lower rate 
of star formation than Fornax or Leo I. Consequently, 
the metal enrichment is less rapid and the SN-induced 
gas blowout is less severe. The resulting MDF has both 
metal-poor and metal-rich stars and is less-peaked than 
fo r the more luminous d Sphs. 

Norris & Bcsscll (1978) first drew attention to the pos- 
sibility that Sc ulptor was chemically inhomogeneous. 
Da Costa (1984) found that the bulk of Sculptor's stars 
are slightly younger than the oldest globular clusters 
(GCs) b ut older than Fornax. W ith HST/WFPC2 pho- 
tometry, iMonkiewicz et al.l (|1999ft found that Sculptor is 
just as old as the GCs. Neither study could determine 
whether the bluer stars were a younger population or 
blue s tragglers from the older population. Mapclli et al. 
(2009) presented evidence that the blue stars are true 
blue stragglers, meaning that Sculptor has only an old 
population. Howev er, old does not neces sarily mean 
single-aged. In fact, iMaiewski et al.l (|1999t ) found that 
Sculptor undoubtedly contains multiple stellar popula- 
tions based on its HB and red giant branch (RGB) mor- 
phologies. The existence of a metallicity spread, the de- 
pression of [a/Fe] with increasing metallicity, and the 
radial change in HB morphology means that star for- 
mation lasted for at least as long as the lifetime of a 
Type la SN and possibly fo r a few Gyr (|Tolstov et al.l 
1200a iBabusiaux et al.112001 . 

Our chemical evolution model conforms to the photo- 
metric description of Sculptor's SFH. According to our 
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Fig. 4. — The observed abundance ratios and the best-fit gas flow and star formation history model for Sculptor. See Fig.[2]for a detailed 
explanation. 



model, Sculptor formed stars for 1.1 Gyr. In fact, one 
of the major advantages of an abundance-derived SFH 
is that it can resolve ages of old populations much more 
finely than a photometrically-derived SFH. As a result, 
we believe our estimate of the star formation duration to 
b e the most precise presen tl y ava ilable for Sculptor. 

lLanfranchi fc Matteuccil ( 20041) also found a chem- 
ical evolution model to match the five stars with 
then- available multi- element abundance measurements 
(jShetrone e t al. 2003). Their model showed a sharp kink 
or knee at the time when Type la SNe ejecta began to 
dilute the [a/Fe] ratio with large amounts of Fe. Our 
model shows a less pronounced knee that occurs at lower 
[Fe/H] and higher [a/Fe] primarily due to our d i fferent 
treatments of the Type la SN DTD. IRevaz et all ()2009l ) 
modeled unpublished abundance measurements by the 
Dwarf Abundances and Radial Velocities Team (DART) 
for Sculptor with a sophisticated hydrodynamical model. 
They found that nearly all of the stars formed between 
10 and 14 Gyr ago, with nearly half of the stars form- 
ing at least 13 Gyr ago. The model supposed that the 
stars formed in about five bursts. It is possible that 
adding burstiness to our model would help to reconcile 
the model with the observed data, such as the peak in the 
MDF at [Fe/H] = -1.3 an d the discrep ancy in [Ca/Fe] at 
high metallicity. However. IRevaz et al.l 's model predicted 
many more stars at [Fe/H] < —3 than we or DART (who 
sample a wider area) observe. A less intense initial burst 
(crudely approximated by the 0.3 Gyr SFR rise time in 
Fig. 01 better matches the low-metallicity MDF. Finally, 



in constructing a c hemical evolution model of Sculptor, 
iFenner et al.1 ()2006| ) found that neutron-capture elements 
contribute significantly to the ability to discriminate be- 
tween different models of star formation. Large, high- 
resolution surveys will add these elements to the dSphs' 
repertoire of abundance measurements. 

Like Fornax, the central regions of Sculptor are dom- 
inated b y a more metal-rich population than the outer 
regions (jBattaglia et al.l [2008T ) . Our sample is centrally 
concentrated in order to maximize the sample size. The 
selection results in a bias toward metal-rich, presumably 
younger stars, possibly shortening the derived the SF du- 
ration compared to what we would deduce from a more 
radially extended sample. 

We also presented Sculptor 's abundance distributions 
in IPaperll (jKirbv et al.1 120091 ) . Minor modifications to 
the abundance measurements (|Paper Hf ) and the restric- 
tion of the plot to points with measurement uncertainties 
less than 0.3 dex in either axis ca use Fig, p] to appear 
slightly different from Figs. 10-12 in lPaper J . The diffcr- 
ences do not affect any of the conclusions of IPaperll . 

3.4. Leo II 

The abundance distributions for Leo II resemble Sculp- 
tor in many ways. The MDF slowly rises to a peak fol- 
lowed by a sharp cut-off, and [a/Fe] declines smoothly 
with increasing [Fe/H]. The best-fit SFH model shows a 
great deal of gas loss, like Sculptor. iBosler et al.l (|2007l ) 
also suggested that Leo II may have experienced more 
intense galactic winds than Leo I due to a lower peak 
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Fig. 5. — The observed abundance ratios and the best-fit gas flow and star formation history model for Leo II. See Fig. [2]for a detailed 
explanation. 



in the MDF. In fact, we find that the mass lost per SN 
(Amt) is higher in Leo II (6.6 x 10 3 M Q SN" 1 ) than in 
Leo I (3.9 x 10 3 M© SN" 1 ). 

Perhaps by virtue of its large Galactocentric distance 
(221 kpc), Leo II h as maintained star fo rmation for 
longer than Sculptor. Mighell fc Rich! ()1996f ) found from 
HST/WFPC2 photometry that the dSph started forming 
stars 14 Gyr ago and continued forming st ars for about 
7 Gyr . In a reanalysis of the same data, lOrban et al.l 
(2008) determined that 30% of Leo IPs stars formed ear- 
lier than 10 Gyr ago and 67% form ed between 5 and 
10 Gyr ago. Shctron e et al.l ([2009) resolved the age- 
metallicity degeneracy in the CMD by using metallicities 
based on spectral synthesis of Kcck/LRIS spectra. They 
found a significant population of stars as young as 3 Gyr. 
However, they pointed out a number of caveats that may 
introduce large errors into their age measurements. 

We derive a star formation duration of 1.6 Gyr. Al- 
though it is the longest duration that we measure for 
the eight dSphs, it does not approach the photometri- 
cally derived durations. The smoothness of the mod- 
eled SFR may mask the true duration of SF. The abun- 
dance distributions — particularly [Si/Fe] and [Ti/Fe]— 
show a smattering of points beyond the main trend line. 
These stars may repres ent stellar populati ons of tempo- 
rally separated bursts. iRevaz et al.l (|2009l ) showed that 
a model with about 13 SF episodes matches the dis- 
persion in [Mg/Fe] at a given [Fe/H] (observations by 
iShetrone et al.|[2009D fairly well. Our model for Leo II, 
like Sculptor, may benefit by adding burstiness. 



3.5. Sextans 

Sextans, Draco, and Ursa Minor form a class of galax- 
ies with similar abundance distributions and SFH mod- 
els. Their MDFs are fairly symmetric (less so for Ursa 
Minor) with a clump of stars at [Fe/H] ~ —3. Their 
[a/Fe] ratios decline smoothly with increasing [Fe/H]. 
The dispersion in [a/Fe] at a given [Fe/H] is fairly small. 
Most of the derived star formation parameters are simi- 
lar (infall normalization, A ln ~ 1.1 — 1.5 x 10 9 Gyr; infall 
timescale, r m ~ 0.2; outflow rate, A out ~ 10 4 Mq SN" 1 ). 

The small bump in the MDF at [Fc/H] ~ —3 deserves 
some discussion because it appears in Sextans, Draco, 
and Ursa Minor. A depression in the MDF appears be- 
tween the bump and the bulk of the MDF. This bump 
might indicate a small, rapid SF burst at early times fol- 
lowed by an epoch of minimal star formation, possibly 
because the SNe from the initial burst blew out the gas. 
When the galaxy reacquired more cool gas, the bulk of 
SF began. The few available [a/Fe] measurements in the 
bump are large, indicating that the stars in the bump 
formed before the onset of Type la SNe. Because our 
model does not permit individual bursts, we can not sup- 
port this speculation beyond our qualitative argument. 

Despit e the low meta ll icity and low luminosity of 
Sextans, Bcllazzin i et al.l ([2001) found that the dSph 
has at least two stellar populations based on its HB 
and RGB morpholo gy. With HST/WFPC2 photometry, 
Orba n et al.l (120081 ) found no stars older than 10 Gyr. 
Lee et al. 120091 measured Sextans 's SFH based on wide 
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Fig. 6. — The observed abundance ratios and the best-fit gas flow and star formation history model for Sextans. See Fig.[2]for a detailed 
explanation. 



field photometry coupled with an algorithm that self- 
consistcntly derives the SFH and chemical evolution of 
the galaxy. They deduced that SF in Sextans occurred 
mainly between 11 and 15 Gyr ago, but some stars 
formed as recently as 8 Gyr ago. However, th ey assumed 
that Sextans is a closed box. In iPaper III we showed 
that the MDF is inconsistent with a closed box. We al- 
low gas to leave the system, which would bring an earlier 
end to SF than in a closed box. As a result, we find a 
SF duration of just 0.8 Gyr. 

3.6. Draco 

Because we conducted a more intense observational 
campaign on Draco than on Sextans, we better sample 
Draco's abundance space. The better sampling does not 
change our qualitative description of the trio comprised 
of Sextans, Draco, and Ursa Minor (see Sec. I3.5|) . The 
metal-rich side of Draco's MDF seems tiered, with fewer 
stars than our model predicts at [Fe/H] = — 1.5 and —1.2. 
The tiers may indicate discontinuous periods of SF. 

As a consequence of its proximity, Draco was one of the 
first dSphs subjected to spectroscopic scrutiny. This sys- 
tem has a stellar mass comparable to globular clusters, 
which are homogeneous in iron-peak elements. There- 
fore, the dis covery of a metal abundance spread within 
this system (iKinman fc Kra ft 1980; Kinman et alJll98ll 
iStetsonl 119841: iSmithl 119841: iLehnert et all 11992ft proved 
to be a notable peculiarity. Furthermore, Draco con- 
tains stars more metal-poor than any globular cluster. 
The first attempt to interpret the metallicity distribu- 



tion within Draco was that of lZinnl (|1978l ). He compared 
metallicities derived for 23 red giants from the Hale 5-m 
multichannel scanner to a chemical evolution model that 
incorporated gas loss (with a rate proportional to the 
SFR) but no gas infl ow. In order to account for the low 
metallicity of Draco, [Zinn (1978) inferred that this sys- 
tem had lost some 90-99% of its initial gas mass. Subse- 
quent spectroscopic and photometric work has more ex- 
tensively documented the MDF and increased the num- 
ber of elements for which abundances have been mea- 
sured (|Shetrone et al.l 119981 120011: lAparicio et al 
i Bellazzini et al.l 120021: iWinnickl 120031: iSmith et al 
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iFaria et al.1120071: IAbiall2008l: iCohen fc Huandl2009h" 

HST/WFPC2 photometry (jGrillmair et all I1998T ) 
and wide- field Isaac Newton Telescope photometry 
(jAparicio et al.l 1200 1| ) showed little evidence for stars 
younger than 10 Gyr in Draco. On the other hand, 
llkuta fc Arimotol ([2002), who also pointed out the simi- 
larities between Sextans, Draco, and Ursa Minor, found 
a lon ger SF duration: betw een 3.9 an d 6.5 G yr. How- 
ever, |Tkuta~^2liniQto|, like iLee et al.l ((2009), assumed 
that a clo sed box w as an adequate description of the 
galaxy. In IPaper III , we determined that failing to ac- 
count for gas outflow overpredicts the peak metallicity 
of the MDF and that failing to account for gas infall 
results in an MDF shape that does not match the obser- 
vations. Our abundance-based SF duration, relaxing the 
closed box assumption, is 0.7 Gyr . Strangely, bas ed on 
the same HS T/WF PC2 data that iGrillmair et~aTl used, 
lOrban et~a l. (2008) determined that half of the stars in 
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Fig. 7. — The observed abundance ratios and the best-fit gas flow 
explanation. 

Draco are younger than 10 Gyr. lOrban et al.l derived 
SFHs for many dSphs, and they did not mention Draco 
explicitly in their text. As a res ult, we do not kn ow why 
t heir SFH diverged from that of iGrillmair et all 

iCohen fc Huand (|2009l ) analyzed high-resolution spec- 
troscopic abundances for eight newly observed stars and 
six stars from the literature. They fit a toy model 
with low- and high-metallicity plateaus in [X/Fe]. The 
low-metallicity plateau has a maximum metallicity of 
[Fe/H] = -2.9 for [Mg/Fe] and -2.4 for [Si/Fe]. We 
do not see a low-metallicity plateau because our sample 
does not include enough metal-poor stars. Instead, we 
observe a smooth, monotonic decline in all four [a/Fe] 
ratios as a function of increasing [Fe/H]. The absence 
of a low-metallicity plateau for the metallicity range of 
our sample suggests that Type la SNe were exploding for 
n early the entire SF lifet i me of Draco. 

IMarcolini etaLl (|2006l I2008D constructed a hydrody- 
namical model of a Draco-like dSph. In order for [a/Fc] 
to drop to 0.2 dex, their modeled dSph must have evolved 
for at least 2 Gyr. However, at small radius — the lo- 
cation of most spectroscopic surveys, including the ma- 
jority of our Draco sample — [a/Fe] does drop to lower 
valu es sooner than i n the dSph as a whole. Nonethe- 
less, IMarcolini et al.l predicted mostly stars with [a/Fe] 
larger than 0.2 dex with a plateau at low metallicity. 
We observe neither of these qualities. Nonetheless, their 
model does qualitatively reproduce important features 
of dSph abundance distributions, including radial gradi- 
ents in both [Fe/H] and [a/Fe], the shape of the MDF, 
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and star formation history model for Draco. See Fig. [2] for a detailed 

and an anti-correlation between metallicity and velocity 
dispersion. 

Finally, we point out that, according to our model, 
Draco lost an enormous amount of gas from SN winds 
during its SF lifetime. lLanfranchi fc Matteuccil (|2007ft 
used Draco and Ursa Minor as case studies in the impor- 
tance of SN winds. One interesting divergence from our 
model is that they found that a wind intensity propor- 
tional to the SFR rather than the SN rate better voided 
the dSph of gas by the present time, in agreement with 
the observed absence of gas. Our different prescription 
for the Type la DTD may mitigate the difference between 
the SFR and SN rate. 

3.7. Canes Venatici I 

Of all of our dSph models, that for Canes Venatici I 
adheres most closely to the observed abundance distribu- 
tions, in part because of the sparse sampling. The MDF 
is a perfect match, and the predicted [a/Fe] line passes 
through the observed locus of points, except for veering 
to slightly high [a/Fe] values at high [Fe/H]. Unfortu- 
nately, only three stars pass the [Mg/Fe] uncertainty cut 
of 0.3 dex. More measurements of [Si/Fe] and [Ca/Fe] 
help us to determine a SF duration of 0.9 Gyr and an 
unusually low SFR exponent of a = 0.36. The weaker 
dependence on gas mass shapes the SFR profile in such 
a way that produces a more symmetric MDF while pre- 
serving a steadily declining [a/Fe] distribution with in- 
creasing [Fe/H]. 

Because Canes Venatici I was discovered recently 
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Fig. 8. — The observed abundance ratios and the best-fit gas flow and star formation history model for Canes Venatici I. See Fig. [2] for 
a detailed explanation. 



( Zuck er et all 20061). few photometric studies exist. 
IMartin et all (|2008bl ) found that the dSph contains 
mostly stars older than 10 Gyr, but 5% of the stars 
could be as young as 1.4 Gyr. iKuehn et al.l ((2008), with 
a shallower CMD, found possible evidence for a pop- 
ulation as young as 0.6 Gyr. They also found three 
candidate anomalous Cepheid variables, indicating an 
intermediate- age population. Because the young popula- 
tion is much smaller than the old population, our chemi- 
cal evolution model and its SF duration should be viewed 
as applicable to the dominant old population. 

3.8. Ursa Minor 

The low-mass Ursa Minor dSph has sometimes been 
studied in comparison with the Draco dSph, in regard 
to bo th its metallicity inhomogeneity and stellar popula- 
tion (lZinnlll98ll: lStetsonlll984t lBelllll985t iShetrorie et al.l 
l200lHBellazzini et al.1 120021 : IWinnickl 120031 : I Abial 120081 ). 
A r elatively small age spread and an ancient mean 
age ([Olszewski fe Aaron^onlll985l : lMighell fe Burkell 19991 : 
iCarrera et al.l 12002ft also makes it also an interesting 
contrast to halo globular clusters. However, spec- 
troscopy has shown that Ursa Minor has a heavy ele- 
ment abundanc e spre a d of more than 1 dex (Zinn 1981; 
Shetrone et all 120011: IWinnickl 120031 : iSadakane et"all 
20041: IGohen fc Huangfe OlO'l even though its stellar mass 



is similar to that of a GC. 

ICudworth. Olszewski, fe Schommen ()1986[ ) conducted 
a photometric survey of Ursa Minor down to the HB. 
With ~ 450 members, they found that the stellar popula- 



tion resembles that of an old, metal-poor GC with a steep 
RGB and a blue horizontal branch The HST/WFPC2 
imaging study of iMighell fe Burke) |l999) confirmed this 
SFH: a single major burst of star formation about 14 Gyr 
ago with a duration of less than 2 Gyr. Our best-fit 
model agrees with these earlier results. From our ob- 
served abundance distributions, we deduce that almost 
all of the star formation in Ursa Minor occurred over an 
interv al of only 0.4 Gyr. In contrast, llkuta fe Arimotol 
(2002) derived an extended period of star formation last- 
ing for ab out 5 Gyr from their closed-box analysis of the 
CMD. In iPaper III we showed that Ursa Minor's MDF 
is inconsistent with a closed box. iCohen fe HuansJ (|2010l ) 
used metallicities from moderate resolution spectra com- 
bined with ages from isochrones to reaffirm that most of 
the stars in Ursa Minor are quite old. 

MDF s have been generated from photometric sur- 
veys by iBellazzini et al.l (120021) and f rom m oderate res- 
olution spectro s copy by IWinnickl (|2003l ). That of 
IBellazzini et al.l (|2002l ) is a good match to our observed 
MDF given in Fig. 9. Both show a sharp rise to a peak 
metallicity of about —2 dcx with a more gradual decline 
towards higher [Fe/H]. The best fit chemical evolution 
model for Ursa Minor produces an MDF that fails to 
match the rapid ri s e seen at [Fe/H] < —2.3 dex. 

ICohen fe Huangl (j2010f ) provided detailed abundance 
analyses for a sampl e of 16 RGB stars, 6 of which ca me 
from earlier wo rk by Shetrone, Cote fc SargenH d200l or 
from ISadakane et al.l (|2004f). Their trends for [Mg/Fe], 
[Si/Fe], [Ca/Fe], and [Ti/Fe] agree qualitatively with 
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Fig. 9. — The observed abundance ratios and the best-fit gas flow 
detailed explanation. 

those found here, but their sample has better coverage 
of the regime [Fe/H] < —2.5 dex, where they found a 
plateau in [a/Fe]. At very low metallicity, [a/Fe] in our 
models reaches highly supersolar ratios, which are larger 
than those observed at the metal-poo r end of the Ursa 
Minor population bv lCohen fc Huaiiel . 

Previous c hemical evolution models of Urs a Minor in- 
clude those of lLanfranchi fc Matteuccl (|2004f ) , who found 
that Ursa Minor has the shortest duration of star for- 
mation of any of the six dSph satellites they studied. 
They deduced that Ursa Minor experienced only a sin- 
gle burst lasting perhaps 3 Gyr, a moderately high 
star formation efficiency, and an intermediate wind ef- 
ficiency In our model the wind efficiency, A out , is the 
highest of all the dSph s in our sample (see Tab. [5]). 
ILanfranchi fc Matteuccf s predicted MDF fails at low 
[Fe/H] as does ours, by being t oo ex tended. In a later 
paper. ILanfranchi fc Matteuccil (|2007fl studied the effect 
of galactic winds. They concluded that a strong galactic 
wind is necessary to reproduce the rather low [Fe/H] of 
the peak of the Ursa Minor MDF, but they still failed to 
reproduce the sudden scarcity of stars more metal-poor 
than t he MDF peak. 

Bot h lMartmez-Delgado et al.l ()2001[ ) and lMunoz et al.l 
(2005) have discovered t idal debris around Ursa Minor. 
As we discussed in Sec. 12.81 (item 9), our observations 
are centrally concentrated and therefore biased toward 
the relatively younger, more metal-rich population that 
is still bound to the dSph. A truly complete analysis of 
Ursa Minor's SFH must also include the tidally stripped, 



<T 100 




<5 ' : 

" -1.0 r 

-1.5 . , 

0.0 0.5 1.0 1.5 

t (Gyr) 



and star formation history model for Ursa Minor. See Fig. [2] for a 

unbound stars. 

4. FURTHER EXPLORATION OF THE CHEMICAL 
EVOLUTION MODEL 

In this section, we explore the parameters of the chem- 
ical evolution model that were previously not allowed 
to vary. Namely, we examine the dependence of the out- 
come of the model on the Type la SN delay time distribu- 
tion, the hypernova fraction, and the metal enhancement 
of supernova winds. We have chosen Sculptor as a case 
study. In each of the following three sections, we alter 
one aspect of the chemical evolution model for Sculptor. 
Then, we use Powell's method to find the combination 
of the six free parameters that maximizes the likelihood 
estimator, as before. A Monte Carlo Markov Chain of 
at least 10 4 trials provides the two-sided 68.3% confi- 
dence intervals for the first two altered models. Table 0] 
compares the results of the new models with the original 
model. 

4.1. Type la Delay Time Distribution 

We have adopted the Type la DTD of iMaoz et~al"1 
(|2C5T0t) . The model is very sensitive to the delay time 
of the first Type la SN to explode after the onset of 
star formation. Unfortunately, this quantity is poorly 
measured. We have chosen 0.1 G yr because that is the 
maximum value that IMaoz et alj 's DTD seems to allow. 
However, the DTD was measured in a range of galax- 
ies with widely varying star formation environments. 
The details of the combined DTD (Fig. [IJ may not be 
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TABLE 4 

Sensitivity of Sculptor Model Parameters to Assumptions 



Parameter Baseline min(t c i c i ay ) = 0.1 Gyr eHN = 0.5 Metal-Enhanced Wind 



A* (10 6 M Q Gyr" 1 ) 0.47±g"?| 0.12±g-g| 0.68±° ^° 0.85 

0-83±oioa °' 93 ±ao9 0.841°$ 0-93 

An (10 9 M Gyr" 1 ) 0.70±°-J* 0.07ig-°J 1.02±g*i« 0.29 

r in (Gyr) 0.27±$™ 0.84+gij* 0.21+g^ 0.13 

Aout (10 3 M a STV" 1 ) 5.36±g-^f 5.38±g-Jf 5.14±g-$ 0.53 

M gas (0) (10 6 M ) 0.50±°i 2 5 O.GOtlH 0.0Ot°:H 0.33 

SF duration (Gyr) 1.05 3.66 0.82 1.26 



Note. — Wc were unable to compute uncertainties for the Metal-Enhanced Wind model because the 
model is numerically unstable to small perturbations. The SF duration is a derived value, not a free 
parameter, and we did not calculate its uncertainty. 
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Fig. 10. — The observed abundance ratios and the best-fit gas flow and star formation history model for S culptor. The dark red (left) or 
heavier (right) lines show the model with the longer minimum Type la SN delay time of 0.3 Gyr (Sec. 1470 . The light blue (left) or faded 
(right) lines show the original value of 0.1 Gyr, as in Fig. [4] 
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single-degenerate Type la SNe will be inhibited at low 
metallicity ([Fe/H] < —1). Nonetheless, the decline 
of [a/Fe] with increasing [Fe/H] in Figs. [3H9] demands 
that some kind of Type la SN explode. Thus, the low- 
metallicity Type la SNe in dSphs may be mergers of 
double-degenerate binaries only. The removal of the 
single-degenerate channel could affect the DTD. 

In order to explore the impact of changing the DTD on 
the chemical evolution model, we have recomputed the 
most-likely model parameters for Sculptor with a mini- 



mum Type la delay time of 0.3 Gyr instead of 0.1 Gyr. 
We did not change the DTD normalization. Figure [10] 
shows the result compared to the original model (Fig. [4]) . 
The abundance distribution is identical except for the 
low- metallicity [a/Fe] plateau, which is flatter for the 
longer delay time because the mass dependence of the 
Type II SN yields is muted. However, the right panel of 
Fig. [10] shows that the SFH has changed dramatically. 
In particular, the timescale of SF has been expanded. In 
fact, the differences in the SFHs can be explained by mul- 
tiplying the time variable in the original model by about 
3.5. The result is less intense star formation over a longer 
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time. In the end, just as many stars are formed and just 
as much gas is blown out as in the original model. 

We conclude that the Type la SN DTD is a major 
uncertainty in our model. The abundance data alone 
does not help to determine the minimum delay time. The 
timescales in our models can be multiplied by a factor 
constrained only by the poorly known minimum Type la 
SN delay time. 

4.2. Hypernova Fraction 

SN 1998bw was immediately identified to be unusual 
because of its association with a gamma ray burst and 
a light curve that sugg ested relativisticall y exp anding 
gas (|Galama et all 119981 ) . Ilwamoto et"all ((19981 ) deter- 
mined that the explosion energy for SN 1998bw was 
about 30 times larger than the average SN. The en- 
ergy of the explosion has con sequences for the nucleosyn- 
thesis. iNomoto et all ([2006) calculated nucleosynthctic 
yields for SNe at a variety of explosion energies. 

One of the fixed parameters in our model is the fraction 
of stars that explode as very energetic hypernovae (chn)- 
We initially chose crn = (no HNe) because it seemed 
to better match the abundance patterns at the lowest 
metallicities (e.g., [Ca/Fe] in Sculptor). In order to ex- 
plore the effect of HNe, we have also found the most likely 
mod el for Sculptor with £hn = 0.5. This is the value 
that INomoto et all chose for their o wn chemical ev o lution 
model of the solar neighborhood. iRomano et all (|2010f ) 
further explored the effect of changing ehn- 

Figure Qj] compares the result of the model with 
£hn = 0.5 with the original model (ehn = 0). The 
abundance distributions are nearly identical except at 
[Fc/H] < —2.3. The model with larger 6hn reaches 
higher [Fe/H] before Type la SNe turn-on. This en- 
sures that the lowest metallicity stars are not polluted 
by Type la SNe ejecta. The result is a plateau in [a/Fe] 
at low [Fe/H]. We further discuss the presence of such a 
plateau in the [Ca/Fe] ratio of Sculptor and the absence 
of plateaus in other dSphs in Sec. 15.1.11 

The effect on the SFH is more noticeable than on the 
abundance distributions. The total star formation dura- 
tion shortens to 0.82 Gyr from 1.1 Gyr. The HN model 
also requires no initial gas, though the original model for 
Sculptor already did not require very much gas. Less gas 
is lost to supernova winds in the HN model 

In conclusion, the inclusion of HNe has a minor effect 
on the abundance distributions and SFH. The most no- 
table result is that very metal-poor stars ([Fe/H] < —2.3) 
in the HN model have [a/Fe] ratios that are inconsistent 
with any amount Type la SN ejecta. Instead, these stars 
incorporate the ejecta of only Type II SNe or HNe. 

4.3. Metal-Enhanced Supernova Winds 

The SNe in our model expel gas without regard to 
its composition. However, SN winds might be ex- 
pected to be more metal-rich than the average gas-phase 
metallicity because metals are more opaque (and there- 
fore more susceptible to radiation pressure) than hy- 
drogren and helium and because the sa me SNe that 
create the metals could b low them away (|Vaderl [1986: 
iMac Low fc Ferraralll999T) . In this section, we explore 
the effect of a metal-enhanced SN wind. We refer 
the reader to iRobertson et all (|2005f) for a more thor- 



ough discussion of a model that included metal-enhanced 
winds from dwarf galaxies. 

We paramaterize the metallicity dependence of the 
wind by fz, which can vary between and 1. Thus, 
we replace Eq. [15] with 
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If fz = 0, then the wind is unenhanced. If fz = 1, then 
the winds expel only metals and no hydrogen or helium. 
For this experiment, we fix fz at 0.01. Although that 
value seems small, the effect on the SFH is dramatic. 

The modeled metallicity distribution (Fig. [12]) does not 
fit the observed distribution as well as for the original 
model. Instead, there is an overabundance of metal-rich 
stars. The metal-rich discrepancy could be mitigated by 
increasing A out (the total amount of gas lost per SN) at 
the cost of worsening the match at intermediate metal- 
licities. The predicted [a/Fe] distributions change only 
at [Fe/H] > —1.2. Metal-enhanced gas loss causes the 
hook back toward lower [Fe/H] in the [a/Fe] diagrams. 
Because the SFR is very low by the time [Fe/H] begins 
to decrease, very few stars are formed during this time. 

The most dramatic effect on the SFH is that much less 
gas is lost over the lifetime of SF in the metal-enhanced 
wind model than in the original model. With an unen- 
hanced wind, Sculptor ejects 1.8 x 10 s M Q of the gas that 
it starts with or accretes. With a metal-enhanced wind, 
that number decreases to 4.5 x 10 6 M . In both models, 
Sculptor forms about 1.2 x 10 6 M Q of stars. The implica- 
tions for galaxy evolution are dramatic. In the first case, 
over 10 8 Mq of gas is required to catalyze star forma- 
tion in Sculptor. Nearly all of this gas is returned to the 
ISM. In the metal-enhanced wind case, star formation in 
Sculptor requires a gas mass of only a few times its final 
stellar mass. The mass of metals returned to the inter- 
galactic medium in both cases is the same, but in the 
metal-enhanced wind model, the metals in the ejected 
gas are much more concentrated. Changes to other as- 
pects of the SFH are subtle. 

We conclude that the amount of metal enhancement 
in the SN blowout dramatically affects the gas dynam- 
ics of the dSph. Even a 1% metal enhancement re- 
duces the total amount of gas required for star forma- 
tion by a factor of 40. However, a model with fz = 
0.01 results in a worse match to the observed metal- 
licity distribution than the original model with an un- 
enhanced wind. A lower, non-zero value of fz might 
produce better agreement with the observed abundance 
data while reducing the amount of gas infall required 
from the unenhanced wind scenario. The literature on 
galactic chemical evolution contains a diversity of SN 
feedback treatments. We refer the reader to the arti- 
cles we have already mentioned (e.g.. iRecchi et aH 



2001 



Lanfranchi fc Matteuccl 120041; IRobertson et all 12005 
Romano et all 2006: M arcolini et alll2008l ) for more thor 
ough treatments. 

5. TRENDS WITH GALAXY PROPERTIES 

We now discuss trends of the abundance distributions 
and derived SFH parameters with observed galaxy prop- 
erties, such as luminosity, velocity dispersion, half-light 



20 



Kirby et al. 




-3.5 -3.0 -2.5 -2.0 -1.5 -1.0 -0.5 0.0 0.0 0.5 1.0 1.5 

[Fe/H] t (Gyr) 



Fig. 11. — The observed abundance ratios and the best-fit gas flow and star for mati on history model for Sculptor. The dark red {left) or 
heavier {right) lines show the model with a hypernova fraction of ejjN = 0-5 (Sec. 14. 21 . The light blue {left) or faded {right) lines show the 
original value of eHN = 0, as in Fig. [4] 



radius, and Galactocentric distance. We show that lumi- 
nosity is the only galaxy property that shows any con- 
vincing correlation with the properties of the abundance 
distributions. 

5.1. General [a/FeJ Trends 

Figure [T3l shows the trend lines of the different element 
ratios with [Fe/H]. The trend line is defined by the aver- 
age of the element ratio, weighted by the inverse square 
of the measurement uncertainties, in a moving window 
of 0.5 dex in [Fe/H]. The moving averages relax the un- 
certainty cut of 0.3 dex (used for Figs. [2HH]) to 1 dex, 
meaning t hat all of the measurements from the catalog 
(jPaper IS ) are included. The bottom panel shows the av- 
erage of four element ratios, which is called ([a/Fe]). The 
weight of the line fades as fewer stars contribute to the 
average near the ends of the MDF. The figure legend lists 
the dSphs in order of decreasing luminosity. For compar- 
ison, some panels of the figure also display the trends for 
the Milky Way ha lo and disk for available element ratios 
(|Venn et alJI2004lfl . 

Fig. [T3] presents the broad trends of the evolution of 

7 The data from I Venn et al.1 12004 | ) is a compilatio n 
of data from the following sources : iBensbv et alj 12003) . 
I Burris et all J200CI'). |E dvardsso n et al.1 111993). IFulbrigM H2000I. 
120021) IGratton fe Snedenl 11983. 11 991. 1994). Hanso netaLl 
11998 | ) , Ilvans et al.1 02003), I Johnson! 120021). I McWillia meTaT] 
|T99a). IM cWilliaml 1199811. [N isscii & Schusteg 1199711 . 
IProchaska et al.l 1200ol). IReddv et all 12003 ). IRvan etall uMB ), 
and Stephens & Bocsgaard (2002). 



[a/Fe] with increasing [Fe/H]. It does not convey the 
width of the dispersion of the [a/Fe] distributions at a 
given metallicity, nor does it show the details at the mar- 
gins of the MDF. The extremely metal-poor stars, which 
represent some of the oldest known stars, are not shown 
in Fig. [H 

5.1.1. Universal Abundance Pattern in dSphs 

The figure does show that the abundance distribu- 
tions of dSphs evolve remarkably similarly. Although 
the dSphs span different ranges of [Fe/H], ([a/Fe]) fol- 
lows roughly the same trend line. This similarity con- 
tradicts the reasonable expectation that different dSphs 
shoul d show a knee in [a/Fe] at different values of [Fe/H] 
("e.g-lMatteucci fc Broc "atol ll99rjtlGilmore fc Wvsel ll991l; 
iTolstov et al.ll2009fl . In fact. iTolstov et al.l didindeed 
find a knee in at [Fc/H] = —1.8 in DART's preliminary 
measurements for [Ca/Fe] in Sculptor. Our measure- 
ments of [Ca/Fe] in Sculptor also show a knee at the 
same metallicity and the same [Ca/Fe]. Ursa Minor pos- 
sibly has a knee in [Ca/Fe] , but with a l ower [Ca/Fe] 
plateau. In agreement with ITolstov et alj 's and others' 
predictions for lower mass systems to experience less in- 
tense SF, Ursa Minor's possible knee occurs at lower 
[Fe/H] than Sculptor's knee. However, the knee is ap- 
parent only in [Ca/Fe] and only in Sculptor and possibly 
Ursa Minor. The element ratios that would better iden- 
tify the onset of Type la SNe, [Mg/Fe] and [Si/Fe], do 
not show a knee for any dSph. 
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Fig. 12. — The observed abundance ratios and the best-fit gas flow an d sta r formation history model for Sculptor. The dark red (left) or 
heavier (right) lines show the model with a metal-enhanced wind (Sec. [4731) ■ The light blue (left) or faded (right) lines show the original 
model with an unenhanced wind, as in Fig. [4] 



The lack of knees for [Fe/H] > —2.5 and the lack of 
low-metallicity plateaus in the [a/Fe] distributions im- 
plies that Type la SNe exploded throughout almost all 
of the SFHs of all dSphs. Of course, the very first 
stars, which have yet to be found, must be free of all 
SN ejecta. The stars to form immediately after the first 
SNe must incorporate only Type II SN ejecta. The very 
lowest metallicity stars in dSphs likely represent this pop- 
ulation. Stars with [Fe/H] > —2.5 formed after the 
Type la SN-induced depression of [a/Fe]. We have al- 
ready explored the possibility of low-metallicity plateaus 
in [Ca/Fe], but we discount the absence of Type la SN 
products as the cause because [Ca/Fe] is the only el- 
ement ratio to show the plateau. We speculate instead 
that metallicity-dependent Type la nuc leosynthesis (e.g., 
iTimmes et alll2003t IHowell et a!1l2OO90 might shape the 
[Ca/Fe] distribution differently from the other element 
ratios. 

High-metallicity plateaus can form when the SF 
achieves a constant rate for a duration long enough for 
the ratio between Types II and la SNe to be constant. 
The SFR would achieve an equilibrium between the pro- 
duction of a elements and Fe. The value of [a/Fe] at 
the plateau depends on the IMF and SN delay time dis- 
tribution. The SF R need not be strictly constant. As 
IRevaz et alj (|2009Q pointed out, a bursty SF profile with 
a high duty cycle can mimic a constant SFR. In that 
case, we would expect a scatter about the mean value 
of [a/Fe] at a given [Fe/H], but the mean value would 



not necessarily evolve with increasing [Fe/H]. We do ob- 
serve high-metallicity plateaus, seen in Fig. [TTJ] The 
trends for [Mg/Fe] and [Si/Fe] do not completely flat- 
ten, but the slopes at [Fe/H] > — 1 are less than the 
slopes at [Fe/H] < -1.5. The trends for [Ca/Fe] and 
[Ti/Fe] do completely flatten for some dSphs. Only 
the more luminous dSphs, which reached metallicities 
of [Fe/H] > —1.2, achieved the high-metallicity plateau. 
The [a/Fe] ratios of Sextans, Draco, Canes Venatici I, 
and Ursa Minor do not flatten. We conclude that dSphs 
with high enough SFRs to reach stellar masses of at least 
10 6 Mq experienced roughly constant SF at late times, 
corresponding to metallicities [Fc/H] > —1.2. 

Beneath the apparently universal path in [a/Fe] -[Fe/H] 
space, the abundance trends vaguely group by luminos- 
ity. Higher luminosity dSphs tend to have slightly higher 
values of [a/Fe] at a given [Fe/H] than lower luminosity 
dSphs. The tracks for Sextans, Draco, and Canes Ve- 
natici I tend to lie below the other dSphs. Fornax and 
Leo I tend to lie above Sculptor and Leo II. These di- 
visions ar e reminiscent of the groupings we proposed in 
iPaper Illl based on MDF shapes. We classified Fornax, 
Leo I, and Leo II as "infall-dominatcd" and Sextans, 
Draco, Canes Venatici I, and Ursa Minor as "outflow- 
dominated." Sculptor sat in its own class. The similar 
groupings based on MDF and [a/Fe] unsurprisingly reaf- 
firm that the SFH shapes both the MDF and the element 
ratio distributions. 

The MW satellite galaxies more luminous than For- 
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Fig. 13. — The moving averages, inversely weighted by measure- 
ment uncertaint y, of abundance r atios for the eight dSphs and for 
the Milky Way (Venn ct al. 2004, who compiled data from the ref- 
erences given in footnote 7). The bottom panel shows ([a/Fe]), 
the average of the top four panels. The line weight is propor- 
tional to the number of stars contributing to the average. The 
legend lists the dSphs in decreasing order of luminosity. Except 
for [Ca/Fe] in Sculptor, the abundance ratios do not show a low- 
metallicity plateau, which indicates that Type la SNe explode for 
nearly the entire duration of star formation. Our data are sparse 
at [Fe/H] < —2.5, and Type la SNe need not explode at times 
corresponding to those low metallicities. Only the galaxies lumi- 
nous enough to reach [Fe/H] > —1 eventually achieve an equilib- 
rium between Types II and la SNe and therefore a plateau at high 
mctallicity. 

nax sample a regime of greater integrat e d star formation 
and higher metallicity. Pompe ia et al.l ([2008D measured 
[a/Fe] for individual red giants in the disk of the Large 
Magellanic Cloud (LMC), and iMucciarelli et al.l (|2008f) 
measured the same for red giants in LMC globular clus- 
ters. The stars span the range —1.2 < [Fe/H] < —0.3 
with one additional star at [Fe/H] = -1.7. The [Ca/Fe] 
ratios of the disk stars decline slightly with increasing 
[Fe/H], but the other element ratios are nearly flat. In 
fact, the LMC stars seem to follow the same [a/Fe] trends 
as Fornax or Leo I, albeit shifted to higher [Fe/H], ex- 
cept for [Ti/Fe]. The average [Ti/Fe] in the LMC is about 



0.1 dex higher than Leo I and 0.3 dex higher than Fornax. 
The Sagittarius dSph also shows a high er average [Ti/Fe] 
than Fornax or Leo I (jChou et al.ll2010f) . Also, [Ti/Fe] in 
Sagittariu s declines wit h increasing [Fe/H] over the entire 
range that lChou et all sampled (-1.5 < [Fc/H] < +0.1). 

The available evidence indicates that the evolution of 
[a/Fe] with [Fc/H] is nearly universal in MW satellite 
galaxies except for [Ti/Fe] at [Fe/H] > -1.3. The av- 
erage values of [Ti/Fe] for the dSphs and the LMC at 
these metallicities vary from about —0.3 (Sculptor) to 
0.0 (LMC and Sagittarius), and the slopes vary from 
A [Ti/Fe] /A [Fe/H] « -0.8 (Sagittarius) to 0.0 (Fornax). 
Ti is both an a element and an iron-group element, and 
it has an appreciable yield fro m both Types II and la 
SNe (|Wooslev & Weaver! [19951) . Therefore, [Ti/Fe] re- 
sponds to changes in the SFR and the IMF differently 
from the "purer" a elements, like Mg and Si. Unfortu- 
nately, our chemical evolution model failed to reproduce 
realistic values of [Ti/Fe] because the theoretical Type II 
SN yields of Ti were too small. We suggest that future 
work explore ratios such as [Mg/Ti] to better understand 
why [Ti/Fe] behaves differently in different dwarf galax- 
ies at high [Fe/H]. 

5.1.2. [Mg/Fe] 

Our data set for the first time has enabled the ex- 
ploration of the bulk properties of [a/Fe] in dSphs that 
span two orders of magnitude in luminosity. In particu- 
lar, Fig. [TBI shows that [Mg/Fe] values higher than in the 
MW are no t unique to the extre mely metal-poor stars in 
dSphs fe.g.. lFrebel et al.ll2010bl ) but also exist in stars of 
more modest metallicity ([Fe/H] < —1.8). 

Factors beyond the SFH may affect the absolute value 
of [Mg/Fe] and other element ratios at low metallicity. 
First, changing the IMF alters [a/Fe] because Type II 
SN yields depend on the mass of the exploding star. 
Second, the early gas mass of the dSph might change 
the shape of the low- metallicity [a/Fe] distribution also 
because SN yields depend on mass. The first SNe in a 
galaxy can more efficiently enrich a small gas mass than a 
large gas mass. Massive SNe explode before less massive 
SNe, and massive SNe generally produce higher [a/Fe]. 
As a result, [a/Fe] at low metallicity could depend on the 
initial gas mass that was enriched by the first SNe. This 
effect possibly explains the larger [Mg/Fe] in dSphs than 
in the MW. We suggest that the stars at [Fe/H] - -2.5 
in dSphs were enriched by SNe of higher average mass 
than the stars at [Fe/H] ~ -2 in the MW. Finally, the 
shape of the abundance distribution might depend on the 
early gas mass because SN yields also depend on metal- 
licity. In addition to sampling higher mass SNe, stars at a 
given [Fe/H] in a lower mass galaxy sample lower metal- 
licity SNe than stars at the same [Fe/H] in a higher mass 
galaxy. 

5.1.3. Unexplained Details 

Many details in Figure [13] defy obvious explanations. 
For example, the [Ca/Fe] ratio is flatter than the other 
element ratios. Sculptor has a strangely large [Ca/Fe] 
at low [Fe/H]. The [Si/Fe] trend for Fornax is above 
the other dSphs' trends, but the other element ratios 
seem consistent. Similarly, the [Ti/Fe] ratio — and only 
[Ti/Fe] — for Leo I lies above the other dSphs. Ursa Mi- 
nor, despite being the least luminous dSph in the figure, 
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has the second largest [a/Fe] at a given metallicity for 
much of the metallicity range. The slope of [Mg/Fe] flat- 
tens for all of the dSphs at [Fe/H] > —1.2, but the slope 
of [Si/Fe] flattens only for Fornax and Leo II. 

We suggest that future work examine the abundance 
catalog in more detail. For example, element ratios with 
a denominator other than Fe could constrain the IMF. 
The predicted yields of [Mg/Si] decrease from +0.2 for a 
progenitor mass of 18 Mp, to — 0-3 for a progenitor mass 
of 40 Mq (jNomoto et al J 120061 ) . Our data set possesses 
the sample size and precision to address such questions. 

5.2. Trends in Chemical Evolution Model Parameters 

We now invoke the best-fit parameters of the chemical 
evolution model in a more quantitative discussion of the 
correlation between abundance distributions and galaxy 
properties. Figure [14] presents the parameters against 
luminosity, line-of-sight velocity dispersion, half-light ra- 
dius, and Galactocentric distance. In addition to the 
model parameters, the bottom row of the figure shows 
the star formation duration, which is a quantity derived 
from the best-fit model, not a free parameter. 

Luminosity can reasonably be expected to show the 
best correlation with quantities related to SF. Of the 
four abscissas in Fig. [TJ1 L is the only one that could 
be predicted from our simple chemical evolution model. 
Roughly, L is the integral of past SF, modulated by the 
reddening and dimming associated with aging. There- 
fore, it is not surprising that the chemical evolution pa- 
rameters vary with L. Although we have plotted the 
SF parameters against L, L is not necessarily the inde- 
pendent variable. Luminosity is a present-day quantity, 
and the stars did not know the final stellar mass of the 
galaxy while they were forming. The SFH determines 
the present luminosity. 

5.2.1. Star Formation Rate Parameters 

The SFR normalization, A„, is roughly constant at 
~5x 10 5 Mq Gyr -1 for galaxies less luminous than 
Leo I. The value roughly doubles for Leo I and increases 
by an order of magnitude for Fornax. The increase in 
A* is expected because a more luminous galaxy must 
have formed more stars than a less luminous galaxy. If 
the SF timescale does not change much with luminosity, 
then the SFR must. We observe that the SF duration 
changes by a factor of about four across the luminosity 
range. Therefore, we estimate a range of 40 in luminos- 
ity The actual L range is 80, but our simple estimate 
ignored the ages of the stellar population and the other 
model parameters which affect the SFR, such as T m . 

The exponent of the SFR law, a, also varies with L. 
If we assume that SFR is proportional to gas volume 
density, then a may indicate the degree to which the 
gas was concentrated in the center of the galaxy. How- 
ever, we find no correlatio n between a and the cqneen- 
tration of the light profiles ()Irwin fe H atzidimi trioulll995l 
not shown in Fig. [T4|) . Our interpretation of a is purely 
speculative because SF is a complex process affected by 
many external factors, such as an ionizing radiation back- 
ground. These factors b ecome more difficult to pr edict 
for smaller galaxies (e.g.. iGnedin fe Kravtsovll2010[) . 

5.2.2. Gas Infall Parameters 



The intensity of infalling gas (or gas cooling to become 
available for SF) drives the SFR. The parameter A m is 
closely related to ^4* . The dSph cannot maintain a high 
SFR without the addition of new gas. Therefore, a lumi- 
nous galaxy must have had large values of both A* and 
A- ln . Alternatively, a luminous galaxy could have started 
its life with a large reservoir of gas. However, in order to 
prevent too many metal-poor stars from forming early, 
new gas must have been added during the SF lifetime. 
The net result is that A*, Ai n , and M gas (0) are highly 
covariant. 

The most likely timescales for gas infall (or cooling) 
vary from 0.17 to 0.42 Gyr. It may be significant that 
none of the timescales exceeds 0.42 Gyr. We propose 
three conjectures. First, T m may reflect the time the 
dSph requires to accu mulate gas. The central densities 
of dSphs are simila r (jMateol Il99l iGilmore et all 120071 : 
IStrigari et al.ll2008h . Therefore, the similar gravitational 
potentials of the dSphs themselves might enforce simi- 
larly small gas accretion timescales. 

Second, the dSphs' environment may set the Tj n 
timescale. Interestingly, ~ 0.1 Gyr was the timescale 
for the Galaxy's monolithic collap se proposed by 
lEggen. Lvnden-Bell. fe S andage (|1962j ). This collapse 
time corresponds to a period when the gas in the vicin- 
ity of the MW was rapidly coalescing into individual 
structures, such as the proto-Galaxy and the dSphs. Af- 
ter 0.1 Gyr, gas accretion would have declined consid- 
erably because the MW and its satellites would by then 
have accreted the bulk of the surrounding gas. In the 
ACDM paradigm, the formation time for a dSph-sized 
dark matter halo is only 0.4 Gyr after the Big Bang 
(jWechsler et al.l 12002ft . Therefore, our most likely gas 
accretion timescales are consistent with both cosmogo- 
nies. 

Third, the time from the formation of the first 
stars to cosmological r eionization is roughly 0.5 Gyr. 
iRicotti fe Gnedinr f2005) referred to all eight of our dSphs 
as "true" or "polluted fossils," meaning that all or most 
of their stars formed before reionization. Our models are 
sensitive to the bulk of the population, and not the few 
younger stars present in most dSphs. Therefore, the best- 
fit values of T m may be probing the pre-reionization SF 
timescale. Fornax must be a exception because the bulk 
of its population formed after reionization. The major- 
ity stellar populations in other dSphs may be fossils with 
SF timescales on the order of the reionization time. The 
(small) dispersion among our r; n values may be a result 
of temp orally protracted, spatially in homogeneous reion- 
ization (|Miralda-Escude et al.l l2000h . However, we note 
that our derived SF durations are longer than 0.5 Gyr 
except for Ursa Minor. To the extent that these dura- 
tions are accurate, we surmise that reionization is one of 
several mechanisms that inhibited SF in dSphs. 

5.2.3. Supernova Winds 

The role of SN feedback f or dSp hs has been emphasized 
repeatedly. iDekel fe Sill (|1986f ) posited that SN feed- 
back regulates the SFR for dwarf galaxies. It can cause 
a terminal wind, or it can blow out gas that is later re- 
accreted. For the smallest galaxies, including the dSphs 
prese nted here, radiation feedback also plays a significant 
role (IDekel fe Wooll2003[ ). The best-fit SN wind intensi- 
ties, A ou %, also show a strong correlation with L. More 
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Fig. 14. — The t rends of the best-fit chemical evolution models with galaxy properties: luminosity (Irwin & Hatzidimitriou 1995; 
[Martinet al. 2008ai), line-of-sight velocity dispersion, two-dimensional projected half-light radius (both from Wolf e t al, 2010|, and ref- 
erences therein), and Galactocentric dis tance. From top to bottom, the parameters are the SFR n orm alization, SFR exponent (Eq. [5J, gas 
infall rate, gas infall timescale fEq. |14> . gas expelled per SN (Eg. 1151 1. and initial gas mass (Eg. 1161 ■ Table [2] gives the same data. The 
bottom row also shows the duration of star formation (see Table [3J, which is a quantity derived from the model, not a free parameter. The 
only galaxy property to show a trend is luminosity. 



luminous dSplis experienced more intense winds. This 
trend is a direct res ult of the m etallicity-luminosity rela- 
tion for d Sphs fe.g. JPaper IIJ l. For reasons discussed in 
iPaper IlH . more intense gas outflow lowers the effective 
metal yield. Therefore, the less luminous, more metal- 
poor dSphs naturally show more gas outflow. However, 
we expected that A out also correlate with the velocity 
dispersion, a measure of the depth of the potential well. 
No such correlation exists. The lack of correlation is puz- 



zling, but the gas blowout depends on the unmeasurablc 
mass density profile at the time of SF and on the loca- 
tions of the SNe within the gravitational potential. 

5.2.4. Galaxy Properties Other Than Luminosity 

The model parameters are insensitive to galaxy prop- 
erties other than L. The velocity dispersions of dSphs 
do not span nearly as large a range as their luminosities, 
which may partly explain the lack of dependence on <7i os . 
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The half-light radius and luminosity together are related 
to the galaxy's surface brightness and stellar density. It 
does not seem that the SF parameters in our model de- 
pend significantly on these quantities. The timescales, 
Ti n and the SF duration, may depend weakly on Galac- 
tocentric distance. The Pearson linear correlation coeffi- 
cient between n n and -Dec is 0.69. Because Ti n basically 
represents the SF duration (the correlation coefficient be- 
tween Ti n and the SF duration is 0.96), this relation may 
indicate that more distant dSphs survive SF-truncating 
inter actions with the MW l onger than closer dSphs. In 
fact. iSilk. Wvse. & Shields! (|1987| ) suggested host galax- 
ies competed with their satellites for gas accretion. The 
more distant satellites, such as dwarf irregulars, success- 
fully accreted more gas to power present star formation 
than the closer satellites, such as dSphs. Orbital his- 
tory would be a better indicator of past interaction with 
the MW. Orbital parame ters based on prop er motions 
are available for F ornax (|Piatek et all I2007D. Sculpto r 
(iPiatek et al.ll20 06). and Ursa Minor (|Piatek et ai]|2005D . 
Soh n et al.l (120071) also constrained the orbit of Leo I 
based on the shape and dynamics of tidal debris. We 
leave orbital analyses for future work. 

We conclude that luminosity is more directly related 
to a dSph's SFH than dynamical or morphological prop- 
erties. The present luminosity can not drive the past 
star formation, but the luminosity does mirror a single 
parameter which determines the SFH. This conclusion 
is simila r to the fundamen tal line for dwarf galaxies de- 
fined bv lWoo et all (|2008t ). They also found that stellar 
mass (closely related to luminosity) is the best predic- 
tor of other dSph properties. However, stellar mass loss 
by tidal stripping may obfuscate the correlation between 
present stellar mass and past star formation. 

6. SUMMARY AND CONCLUSIONS 

We have made a first attempt at quantitative chemical 
evolution models for the large sample of multi-element 
abundanc e measur ements for MW dSphs that we pub- 
lished in iPaper III . Our simple model is a significant 
improvement to the analytical mode ls of the m etallic- 
ity distributions that we explored in IPaper III We fit 
the MDF and [a/Fe] distribution simultaneously to de- 
rive the SF and gas flow histories of each of eight dSphs 
spanning about two orders of magnitude in luminosity. 
Our model produces reasonable fits to the abundance 
distributions of dSphs whose color-magnitude diagrams 
show that most or all of their stars are older than 10 Gyr. 

We draw the following conclusions from our models 
and from the general trends in abundance distributions 
(Fig. US]): 

1. The [a/Fe] ratios evolve with metallicity along 
nearly the same path for all dSphs. The average 
value of [Mg/Fe], [Si/Fe], [Ca/Fe], and [Ti/Fe] drop 
from +0.4 at [Fe/H] = -2.5 to 0.0 at [Fc/H] 

— 1.2, where the slope flattens. 

2. No low- metallicity plateaus or knees exist in [a/Fe] 
vs. [Fe/H] space for any dSph at [Fe/H] > —2.5. 
We conclude that Type la supernovae contributed 
to chemical evolution for all but the most metal- 
poor stars. 



3. The [Mg/Fe] ratio in dSphs exceeds that of the 
Milky Way at [Fe/H] < -1.8. We suggest that 
the abundance ratios of stars in low-mass systems 
are more sensitive to the mass and metallicity de- 
pendence of Type II supernovae yields than stars at 
the same metallicity in higher-mass systems, such 
as the progenitors of the inner MW halo. 

4. The dSphs may be grouped based on their [a/Fe] 
distributions into roughly the same groups that 
we de fined base d on their metallicity distribu- 
tions (jPaper HID . The more luminous dSphs 
have infall-dominated MDFs and slightly higher 
([a/Fe]) at a given [Fe/H]. The less luminous dSphs 
have outflow-dominated MDFs and slightly lower 
([a/Fe]) at the same [Fe/H]. 

5. Some SF model parameters correlate with present 
luminosity, but not with velocity dispersion, half- 
light radius, or Galactocentric distance except for 
a possible correlation between gas infall timescale 
and D GC - 

6. The gas flow histories for all dSphs except For- 
nax are characterized by large amounts of gas loss, 
probably driven by supernova winds. Less lumi- 
nous dSphs experienced more intense gas loss. 

7. Allowing supernova winds to be metal-enhanced 
drastically reduces the amount of gas infall and 
outflow required to explain the observed abundance 
distributions. 

8. The gas infall timescale does not exceed 0.42 Gyr. 
This possibly reflects the amount of time ancient 
stars had to form before reionization ended star 
formation. 

9. The derived star formation timescales are ex- 
tremely sensitive to the delay time for the first 
Type la SN. Increasing the delay time from 0.1 Gyr 
to 0.3 Gyr results in a star formation duration in 
Sculptor inflated by a factor of 3.5. 

10. The presence of bumps in the MDFs and stars with 
[a/Fe] ratios far from the average trend lines sug- 
gests that the SFHs of dSphs were characterized 
by bursts, which are not included in our model. 
Bursts are a common feature of more sophisticated 
models. 

Some of our conclusions (5-10) depend on the realism 
of our chemical evolution model. Many more sophisti- 
cated models exi st, and w e encourage their application 
to our data set. IPaper l3 contains the complete abun- 
dance catalog. 

The major strength of the present work is that we ap- 
ply the same model to a homogeneous data set of hun- 
dreds of stars in each of eight dSphs. The sample size 
and diversity of galaxies has allowed us to present an 
overview of chemical evolution in dwarf galaxies. We 
have discovered patterns not apparent in previous data 
sets due to small samples or lack of diversity among the 
well-sampled galaxies. In particular, we have shown that 
[a/Fe] distributions of dSphs do not form a sequence of 
knees corresponding to the metallicities at which Type la 
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supernovae began to explode. Instead, the [a/Fe] pat- 
terns of all dSphs are largely the same, but different 
dSphs sample different regions in metallicity. 
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